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ABSTRACT 


In this investigation, a comprehensive Vacuum Assisted Resin Transfer Molding (VARTM) 
process simulation model was developed and verified. The model incorporates resin flow 
through the preform, compaction and relaxation of the preform, and viscosity and cure ki- 
netics of the resin. The computer model can be used to analyze the resin flow details, track 
the thickness change of the preform, predict the total infiltration time and final fiber vol- 
ume fraction of the parts, and determine whether the resin could completely infiltrate and 
uniformly wet out the preform. 

Flow of resin through the preform is modeled as flow through porous media, Darcy’s 
law combined with the continuity equation for an incompressible Newtonian fluid forms the 
basis of the flow model. During the infiltration process, it is well accepted that the total 
pressure is shared by the resin pressure and the pressure supported by the fiber network. 
With the progression of the resin, the net pressure applied to the preform decreases as a 
result of increasing local resin pressure. This leads to the springback of the preform, and is 
called the springback mechanism. On the other side, the lubrication effect of the resin causes 
the rearrangement of the fiber network and an increase in the preform compaction. This 
is called the wetting compaction mechanism. The thickness change of the preform is deter- 
mined by the relative magnitude of the springback and wetting deformation mechanisms. 
In the compaction model, the transverse equilibrium equation is used to calculate the net 
compaction pressure applied to the preform, and the compaction test results are fitted to 
give the compressive constitutive law of the preform. The Finite Element/Control Volume 
(FE/CV) method is adopted to find the flow front location and the fluid pressure. The code 
features the ability of simultaneous integration of 1-D, 2-D and 3-D element types in a single 
simulation, and thus enables efficient modeling of the flow in complex mold geometries. 

VARTM of two flat composite panels was conducted to verify the simulation model. The 
composite panels were fabricated using the SAERTEX multi-axial warp knit carbon fiber 
fabric and SI-ZG-5A epoxy resin. Panel 1 contained one stack of the carbon fabric, and 
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Panel 2 contained four stacks of the fabric. The parameters verified included the flow front 
location and preform thickness change. For Panel 1, the flow front locations were accurately 
predicted while the predicted resin infiltration was much slower than measured for Panel 
2. The disagreement is attributed to the permeability model used in the simulation, which 
failed to consider the interface flow in the unstitched preform containing more than one stack 
of the fabric under very low compaction force. The predicted transverse displacements agree 
well with the experimental measurement qualitatively, but not quantitatively. The reasons 
for the differences were discussed, and further investigations are recommended to develop a 
more accurate compaction model. 

The simulation code was also used to investigate the VARTM of a new form of sandwich 
structure with through-the-thickness reinforcements, which is being considered for use in 
primary aircraft structure. The infiltration of three foam core sandwich preforms with dif- 
ferent stitch densities was studied. The objective of the study was to determine whether the 
preforms could be completely infiltrated and how the stitch density affects the infiltration 
process. The visualization experiments were conducted to verify the simulation. The model 
accurately predicted the resin infiltration patterns. The calculated filling times underpre- 
dicted experimental times by 4 to 14%. The model revealed the resin flow details and found 
that increasing the stitch spacing shortens the total filling time, but increases the nonuni- 
formity of the flow front shape. Extreme nonuniformity of the flow front shape could result 
in the formation of the voids. 
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Chapter 1 


Introduction 


The Resin Transfer Molding (RTM) process is a cost-effective fabrication method for the 
manufacture of polymer composites. In a traditional RTM process, catalyzed thermosetting 
resin is injected into an enclosed metal mold containing a previously positioned reinforcement 
preform. The preform is compacted to the specified fiber volume fraction when the matched 
metal mold is closed. The resin wets out the fiber until the mold is filled, and the part is 
then cured inside the mold. A schematic diagram of the process is shown in Fig. 1.1. 

RTM offers several advantages over other composite fabrication methods such as auto- 
clave and compression molding of prepreg tape laminates. First, high fiber volume fraction 
composites can be fabricated with low void contents. Second, parts with highly complex 
shapes can be molded by incorporating many components into a single preform. This helps 
to reduce the cost and weight of the structure. Third, hand lay-up of prepreg tape is elim- 
inated. Production rates are increased and operating costs are reduced. Finally, RTM is a 
closed mold process that reduces the workers’ exposure to harmful volatiles, i.e., styrene, as- 
sociated with many of the room temperature processing resins. However, the matched metal 
tooling used in the RTM procedure is expensive, and the tooling design becomes difficult 
when fabricating large and complex shaped parts. 
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Figure 1.1: RTM processing schematic diagram 

Vacuum Assisted Resin Transfer Molding (VARTM) has been developed as a variant of 
the traditional RTM process to reduce the cost and design difficulties associated with large 
metal tools. In VARTM, the upper half of a conventional mold is replaced by a vacuum bag. 
This eliminates the need for making a precise matched metal mold as in the conventional 
RTM process. Fig. 1.2 shows a generalized schematic diagram of the VARTM process. First, 
the vacuum pump is turned on to expel air from the preform assembly. After the system has 
been equilibrated and all air leaks have been eliminated, the resin is allowed to flow into the 
preform. A pressure of 1 atm provides both the driving force for the resin to impregnate the 
reinforcement and the compression force to compact the preform to the desired fiber volume 
fraction. The vacuum is left on until the resin has completely gelled. The part may then 
be cured at room temperature or in an oven. Due to the low injection pressure (1 atm), a 
resin distribution medium with high permeability is often incorporated into the vacuum bag 
lay-up to facilitate the resin flow in the preform. Because the distribution medium provides 
a low resistance resin pathway, the resin easily fills the distribution medium, and then leaks 
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Figure 1.2: VARTM processing schematic diagram 

downward and wets the preform out. The use of distribution medium reduces the processing 
time significantly, and ensures complete wet-out of the preform. 

VARTM offers numerous cost advantages over traditional RTM, such as lower tooling cost 
and shorter start-up time. Because the top half of the tool is a flexible vacuum bag, a matched 
metal mold is no longer required. This makes VARTM attractive for the manufacture of large- 
scale components. Along with the advantages, VARTM also has inevitable disadvantages. 
It is a relatively complex process to perform well. The flexible nature of the vacuum bag 
makes it difficult to control the final thickness of the preform, and thus, the fiber volume 
fraction of the composite. 

Due to the complex nature of the VARTM process, the trial and error method is inefficient 
and expensive for the process design and optimization. The objective of this study was to 
develop and verify a comprehensive VARTM simulation model as a cost effective design tool. 
The model includes three submodels: the viscosity and cure kinetics of the resin, the resin 
flow in the porous reinforcement preform, and the compaction of the preform during the 
infiltration procedure. The modeling software, 3DINFIL 5.0, has been developed to analyze 
the flow details, track the thickness change of the preform, predict the total infiltration time 
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and final fiber volume fraction of the parts, and determine whether the resin completely 
infiltrates and uniformly wets out the preform. The accuracy of the model was verified by 
comparing experimental results from VARTM of flat preforms and sandwich structures with 
model predictions. 

A review of recent research in modeling of the VARTM procedure is presented in Chap- 
ter 2. Chapter 3 details the resin viscosity and kinetics, flow, and preform compaction 
models. The numerical technique used for the flow model, the finite element/control volume 
(FE/CV) method is discussed in Chapter 4. Comparisons between the model predictions 
and experimental results are shown in Chapter 5 and Chapter 6. Conclusions and recom- 
mendations for future work are given in Chapter 7. 
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Chapter 2 


Literature Review 


2.1 Resin Flow 

Complete filling of the mold with adequate wetting of the fibrous media is the primary objec- 
tive in resin transfer molding (RTM) and vacuum assisted resin transfer molding (VARTM). 
Incomplete impregnation of the reinforcement preform results in expensive defective parts, 
which must be scrapped. Therefore, the simulation of the resin flow through the preform 
is an important step in modeling these composites manufacturing processes to ensure that 
the process can be finished before the resin gels. Simulation codes can be used to track 
flow fronts and estimate the fill times. Modeling the resin flow during the impregnation 
processes provides insight into the process physics and highlights potential problems before 
production. In addition, flow prediction enables optimization of the process variables, such 
as vent location, and injection port location and type. 

Many computer models have been developed to simulate the mold filing process in RTM. 
All these studies treat the propagation of the resin through the reinforcement preform as 
an incompressible fluid through a porous media. Darcy’s law is used to describe the flow. 
Darcy’s law in conjunction with the continuity equation forms the basis for the numerical 
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simulation of the infusion process. Different numerical techniques have been used to model 
this flow, such as finite difference method [1, 2], boundary element method [3], and finite 
element method [4-11]. Of these, the Finite Element/Control Volume (FE/CV) method 
has been widely used because of the lower computational cost. In this approach, the mold 
is first discretized into finite elements. A control volume is constructed around each node 
by subdividing the elements into smaller volumes. The pressure distribution in the resin 
filled region is obtained by use of the finite element method and Darcy’s law is applied to 
calculate the flow rates. The concepts of the control volume and fill factor are then employed 
to track the location of the flow front [4], The advantages of FE/CV method is that the flow 
simulation for complex geometries can be done rapidly and efficiently without the need of 
remeshing the resin-filled domain. 

Most RTM processes deal with parts of a shell-like geometry. Because the part thickness 
is much smaller than its other dimensions, resin flow in the thickness direction has been 
neglected and two-dimensional RTM flow models have been developed by several indepen- 
dent research groups [5-7]. The calculated domain is a three-dimensional (3D) geometry, 
but the flow in the thickness direction is neglected, therefore these models are often called 
2|-dimensional (2|D) flow models. The major advantage of the 2|D model is that the 
computational cost can be greatly reduced with acceptable simulation accuracy achieved. 

Bruschke and Advani [6] created a two-dimensional FE/CV model to predict the resin flow 
pattern for mold filling in anisotropic media. A key feature of the simulation is the ability to 
model the influence of non-isothermal permeability of the fiber preform on the filling process. 
The results of the simulation are compared with experiments performed in a flat rectangular 
mold. A Newtonian fluid was injected into the mold at the center with the pressure at the 
inlet held constant. Four sets of experiments were conducted. In the first set of experiments, 
a lay-up of nine unidirectional mats with [0, ±45, 90, 0, 90, ±45, 0] orientations was resin 
infused to study the effect of stacking layers with different orientations. The second set of 
experiments used a round insert to observe the influence of geometric complexity on the flow 
field. In the third set of experiments, the preform consists of two pieces of unidirectional 
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mats put together. The local permeability at the interface between the pieces is different 
from the bulk material and the preform becomes non-homogeneous. The interface between 
the pieces was called the cutline. This set of experiments was designed to investigate the 
influence of local non-homogeneity along the cutline. The last set of experiments examined 
the resin flow in a mold with two different gap-heights. In all the experiments, the upper 
half mold tool plate was made of glass for the purpose of flow visualization. The flow front 
locations were recorded on video tape to compare with the numerical results. In most cases, 
the predicted and experimental flow front corresponded well. However, the flow fronts in the 
experiment moved faster along the mold walls and the cutline. This showed the existence of 
low-resistance flow paths at these locations. This phenomenon is called race-tracking. 

In most RTM process, resins injected into the mold cavity are highly reactive. The 
resin kinetic and viscosity changes during the infusion process are significant and need to 
be considered in the mold filling simulation. A computer simulation package including resin 
flow, heat transfer and chemical reaction based on the FE/CV method was developed at the 
Ohio State University [7,8]. It can be used to simulate a two-dimensional non-isothermal 
molding process. 

Calhoun et al. [9] conducted experiments at the Dow Chemical Company to verify the 
RTM flow model developed at Ohio State University. The experiment design was completely 
decoupled from the model development to eliminate any possible bias. Rectangular panels 
were injected using RTM process. Two different injection schemes were studied. The first 
scheme was designed to form a one-dimensional flow with a line source on one end of the 
panel and a line sink on the other. The second flow scheme used runners around the preform 
edge to produce a converging flow front to a point sink at the center. In the experiments, 
pressure transducers were used to monitor the pressure change during the infusion procedure 
at different locations in the preform. It was found that the mold filling time could be 
predicted within experimental error but the pressure distribution needed to be scaled to 
match the experiment data. The difference between experiment and simulation was likely 
due to the local variations of the preform permeability. It was also demonstrated that the 
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flow modeling was extremely sensitive to the preform porosity and permeability. The race- 
tracking phenomenon was observed in the one-dimensional flow experiment and different 
edge permeabilities were tried to simulate the phenomenon. When the edge permeability 
was set to be 100 times the bulk permeability, the simulated and experimental pressure 
profiles agreed well with each other. 

Two-dimensional flow models work well to simulate the infusion of parts with small 
thickness. However, thicker parts require extensive resin flow through the thickness and thus 
the three-dimension effect in the mold filling process must be considered. Lim and Lee [10] 
developed a computer code to simulate the three-dimensional RTM mold filling process, 
including the temperature and cure effects. The validity of the numerical code was verified 
using simple problems with analytical solutions known. Experiments were also performed 
to further validate the numerical results. In the experiment, a cubic shaped preform was 
injected from the inlet gate at the center of the bottom surface. An optical fiber was used to 
monitor the resin flow inside the thick preform. The simulated resin flow front location as a 
function of time was compared with the experimental data. The agreements were found to 
be satisfactory. Numerical simulations for the filling process of complex geometries, such as 
an automobile headlamp bezel with a membrane and centrifugal pump casing, were done to 
illustrate the effectiveness of the numerical code. 

In order to model the infusion process for a variety of conditions of industry interest, 
such as runner distribution systems, shell geometries, and 3-D flows, Phelan [11] developed 
a computer simulation model which incorporated a variety of element types. The model 
featured the ability of simultaneous integration of 1-D, 2-D, and 3-D element types in a 
single simulation. The model used the Flow Analysis Network (FAN), which is referred as 
FE/CV method by other researchers, to track the movement of the flow front. The model 
predicted the flow front position as a function of time and the pressure distribution during 
the filling process. The filling of a complex automotive part, the front-end structure of a Ford 
Escort was presented. The simulation proved to be a useful tool to quickly screen potential 
gating problem and control the pressure levels during the filling process. 
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Researchers have employed computer simulation models as a design tool to develop new 
injection schemes for the RTM procedure [12,13]. Chan and Morgan [12] proposed a multi- 
ple port resin injection scheme to reduce the mold filling times. The ports were opened in 
sequence during RTM. A one-dimensional flow model was developed to illustrate the advan- 
tages of the proposed method. Results from the numerical study showed that the progressive 
opening of the injection gates reduced the filling time remarkably. Kang et al. [13] investi- 
gated two different multiple gate injection schemes. One is to open the gates in a sequential 
manner as the resin front advances. The other is to open all the gates simultaneously and 
control the pressures of each gate individually. Mold filling was simulated numerically using 
the FE/CV method. Experiments were conducted to verify the simulation results. The resin 
flow front shapes observed from the experiments compared well with the numerical results. 
Numerical simulation for the mold filling of the front panel of a tourist bus was performed 
to illustrate the efficiency of the proposed injection schemes. Both of the controlled multiple 
gate injection methods were proven to significantly reduce the injection time as well as to 
prevent the formation of air bubbles due to air entrapment. 

The Finite Element/Control Volume (FE/CV) method has been successfully used to 
track resin flow in the RTM process. However, this popular numerical scheme also has some 
inherent drawbacks. First, the exact location of the flow front is difficult to determine as a 
consequence of the fixed mesh system. This ambiguity in the resin flow front location limits 
the accuracy of FE/CV method [14], Mass conservation problems have also been reported 
with the use of this numerical approach [11,15]. Researchers have addressed these numerical 
problems and put forward methods to improve the conventional FE/CV method [14,16]. 

Kang and Lee [14] proposed a new algorithm, referred to as the Floating Imaginary 
Nodes and Elements (FINE) method, to get a smoother flow front and reduce the error in 
the pressure at the flow front of a FE/CV simulation of resin flow in a porous medium. 
With the FINE method, imaginary new nodes were added at the estimated flow front and 
the flow front elements were divided into two separate regions: the area of resin and the area 
of air. Thus, the flow front element was refined in an adaptive manner at each time step. 
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Two examples were studied to evaluate the validity of the FINE method. First, a square 
mold with a circular hole at the center was considered. Resin was injected from the center 
of the bottom edge of the square mold at a constant pressure. In the second case, resin was 
injected at constant flow rate into the square mold through a circular hole at the center. 
The simulations by the FINE method were compared with the analytical solutions and the 
numerical results by the conventional FE/CV method. It was shown that the FINE method 
improved the solution accuracy considerably with a marginal increase of the computational 
cost. 

Violation of the law of mass conservation was reported with the use of the FE/CV 
approach. When the Galerkin formulation is used with isoparametric finite elements to 
obtain the pressure field, the balance of resin mass cannot be achieved [11]. A numerical 
model, ‘the flow redefinition model’, was proposed by Joshi et al. [16] to eliminate the 
problem. In the model, the intra-control volume and inter-control volume flow balances 
were enforced by arranging the control volumes in a descending order of their representative 
pressure and adjusting the outflow to match the inflow for all the control volumes in the 
same sequence. To validate the model, the authors studied the central injection mold filling 
problems where analytical solutions exist. Significant improvements in the fill time prediction 
were achieved using the flow redefinition model. 

All the numerical simulations mentioned above use the conventional physics, which is 
based on Darcy’s law and the continuity equation for the incompressible fluid. The flow of 
resin through the fiber packed mold cavity is modeled as the saturated flow of a Newtonian 
liquid through a single scale porous medium. However, it was noted that for the woven, 
stitched, or braided fiber mats, the preform behind the flow front was not fully saturated and 
the conventional physics is inadequate to simulate the resin impregnation process [17-19]. 
Simple, one-dimensional mold filling experiments were conducted with resin injected at a 
constant flow rate. In such experiments, it was found that the inlet pressure history deviated 
significantly from the linear plot predicted by the conventional physics. This anomalous flow 
behavior observed in the preforms made of the woven, stitched, or braided fiber mats was 
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attributed to the dual-scale medium with the presence of large gaps between dense fiber 
tows [17,18]. 

A two layered model [20, 21] was proposed to simulate the unsaturated flow in woven 
fiber preforms. The movement of the resin through such fiber mats was modeled as the flows 
through the inter-tow and intra-tow pore spaces separately, and these flows were coupled 
by a sink term, which conserved the global mass. A dimensionless number, ‘pore volume 
ratio’, was proposed to quantify the strength of the ‘sink’ effect. The pore volume ratio 
was defined as the ratio of pore volumes in the fiber bundles and the gap regions between 
the bundles. The model recreated the non-linear drooping inlet pressure history of the 
woven mats obtained in the experiments. However, the effect of the dimensionless number 
proposed in the model, ‘pore volume ratio’, could not be validated empirically. The study 
of the unsaturated flow in the woven, stitched, or braided fiber mats is still in its infancy, 
and more research efforts are needed to advance this study [19]. 

Vacuum Assisted Resin Transfer Molding (VARTM) is a variant of traditional RTM 
process in which one of the solid tool faces is replaced by a flexible vacuum bag. The 
process is also named as Vacuum Bag Resin Transfer Molding (VBRTM), Vacuum Bag 
Molding (VBM), and Resin Infusion under Flexible Tooling (RIFT), etc. Williams et al. [22] 
reviewed the progress of VARTM from its first development as the Marco method in 1950 
to the Seemann Composites Resin Infusion Manufacture Process (SCRIMP) today. A series 
of SCRIMP patents were awarded in the USA in the early to mid nineties [23-25]. The 
process is simple and resembles other vacuum infusion techniques in that the laminate is 
contained under a vacuum bag made of air impervious films such as nylon or silicone film, 
and resin is drawn in under vacuum. The novel aspect of SCRIMP is the use of patented resin 
distribution systems to facilitate mold filling and improve production efficiency. Generally, 
there are two types of resin distribution systems to facilitate the resin infusion in the SCRIMP 
process. One uses a distribution medium with high permeability placed on top of the fiber 
preform. The distribution medium is also referred to as the high permeability layer (HPL) 
or the liigh-permeable medium (HPM). Resin flows easily through the distribution medium 
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and then wets out the reinforcement preform from above. The use of a high-permeable 
distribution medium creates a transverse (i.e. , top-to-bottom) flow in the fiber mat, which 
reduces the resin infusion time. The other type of SCRIMP uses grooves cut on a low-density 
core to distribute the resin. The preform is placed on top of the core and resin infuses the 
entire structure from the grooves. The core is part of the final product. The permeability 
of the grooves is much higher than that of the distribution medium, thus the part can be 
infiltrated more quickly. However, SCRIMP based on the distribution medium has attracted 
more attention in industrial applications and research investigations because of the ease to 
process and design flexibility. 

Tari et al. [26] investigated the effect of the high permeability layer (i.e., distribution 
medium) on the resin flow field in VARTM through analysis, simulation and experiment. 
In the analysis, it was assumed that the shape of the flow front does not change as it 
travels down the mold. Fig. 2.1 shows a schematic diagram of the flow front associated 
with the use of a high permeability layer (HPL). In the flow front region, the pressure 
through the thickness of HPL was assumed constant and in the fiber mat, the flow in the 
in-plane direction was neglected compared to the transverse flow. For the flow front region, 
application of Darcy’s law and the continuity equation to the flow in the HPL and fiber mat 
yielded a closed form solution for the flow front shape. Assuming no cross flow between 
the HPL and fiber mat in the region behind the flow front, the location of the flow front 
was solved as a function of time. The 2-D finite element model, RTMSIM developed at 
UCLA, was employed to simulate the resin flow in the cross-section parallel to the flow 
advancement. In the experiments, bi-directional fiber mats with different thicknesses were 
resin infiltrated in a rectangular Plexiglas mold. A commercially available loose weave nylon 
cloth was used as the HPL. The side and top views of the flow front were recorded using 
two video cameras. The experiment results validated the assumption that the flow front 
shape does not change during the mold filling procedure. The numerical simulation gave 
accurate results as to the shape of the flow front. The resin flow in the HPL led the flow 
in the fiber mat. The lags observed experimentally were compared to those predicted by 
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Figure 2.1: Schematic diagram of flow front associated with using a HPL 

the theory and simulation. The simulation results were within 13 percent of the measured 
values. The model overpredicted the length of the lag because the theory did not consider 
the contribution to the fiber saturation from in-plane flow. This study showed that the use 
of HPL created a transverse flow in the fiber mat, and thus greatly reduced mold filling time. 

Hsiao et al. [27] proposed a more accurate closed form solution for flow in the distribution 
medium and fiber mat during the VARTM process by taking the in-plane flow in the fiber 
mat into consideration. The flow was divided into a saturated region with no cross flow and 
a flow front region where the resin was infiltrating into the preform from the distribution 
medium. The flow front region moved with a uniform velocity and kept a steady shape. The 
law of conservation of mass and Darcy’s law were applied in each region. It was found that 
the flow front was parabolic in shape and the length of the saturated region was proportional 
to the square root of the time elapsed. The results obtained were compared with full-scale 
3D numerical simulations and an error analysis carried out. The value of the parameter, 
e = K‘) X x h‘> 2 / K-iyy d' 2 ■ was a measure of the error between the analytical model and the 
numerical simulations. In the formulation, K- Ixx and K^yy are the permeabilities of the fiber 
mat in the direction of flow advancement and thickness direction, respectively. /i 2 represents 
the thickness of the fiber mat and d is the flow front length (Fig. 2.1). The error in the 
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total infiltration time remains small as e increased. However, the error in the flow front 
length grows rapidly while e increases. The reason for this is that the theory neglected the 
transition region between the saturated region and the flow front region. The contribution 
of this region to the through-thickness flow was effectively lumped into the flow front region, 
and thus gave rise to the error. As the transition region increases in significance with 
increasing e, it increasingly affects the flow front region and the error in the flow front length 
grows rapidly. A parametric study was conducted to investigate the influence of a number 
of process variables on the fill time and flow front velocity. It was found that the fill time 
decreases and the velocity increases when the thickness of the distribution medium relative 
to that of the preform increases. The porosity of the distribution medium was found not to 
affect the filling time and the flow front velocity significantly. The flow front velocity was 
found to be much more sensitive to the permeability of the distribution medium than that 
of the preform. 

Sun et al. [28] investigated the flow mechanism in the SCRIMP process based on the dis- 
tribution medium through flow visualization experiments and numerical simulations. Flow 
visualization of SCRIMP mold filling was carried out under various molding conditions. 
It was found that the resin flowed faster in the distribution medium than in the fiber rein- 
forcement, and the flow front lead-lag remained nearly constant during the entire mold filling 
process. Experiments were conducted to evaluate the importance of the distribution medium 
in mold filling. Two preforms were injected under similar conditions. The only difference 
was that the distribution medium was placed on top of the preform in one experiment while 
no distribution medium was used in the other experiment. The use of distribution medium 
reduced the filling time by 85 percent. It was very clear that the distribution medium is 
essential to SCRIMP mold filling. A 3D FE/CV model was employed for the numerical 
simulation of the mold filling process. The numerical results for the filling time compared 
well with the experimental measurements. A parametric study was carried out to investigate 
the influence of the flow properties of the distribution medium and preform on the filling 
time. The simulation results indicated that the mold filling time was a strong function of 
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the permeability of the distribution medium, but was not sensitive to that of the preform. 
Experimental results and the 3D FE/CV simulation demonstrated that the in-plane flow 
in the preform could be neglected when the permeability of the preform is much less than 
that of the distribution medium and the length along the flow direction is much larger than 
the thickness of the composite part. Based on this observation, a leakage flow model was 
presented to consider the flow as a 2D in-plane flow in the distribution medium with a sink 
term to account for the resin leakage into the preform. The leakage flow model reduced the 
computational cost by more than 90 percent and yet provided a reasonable estimation of the 
mold filling time in SCRIMP based on the distribution medium. 


2.2 Compaction of the Preform 

Investigation of resin flow in the porous preform in the VARTM process helps industry to 
ensure that complete impregnation of the preform before resin gellation occurs. Another 
major concern of the manufacturers is the prediction of the final thickness and thereby, the 
fiber content of the composite laminate, which is governed by the compaction behavior of 
the preform. Due to the flexible nature of the vacuum bag, there is no direct control over 
the thickness or fiber volume fraction of the composite part in the VARTM procedure. The 
compaction of the reinforcement preform is complex and depends on the compressibility and 
relaxation of the reinforcement under pressure, and the interaction between the reinforce- 
ment and the resin flow. Researchers [29, 30] have conducted experiments to investigate 
the mechanisms governing compaction of the preform during the VARTM process. The 
experimental observations are interesting and elaborated in the following paragraphs. 

The flexible nature of the vacuum bag, coupled to the varying pressure inside the mold 
cavity, results in the laminates with a nonuniform final thickness. Rigas et al. [29] investi- 
gated variations of the thickness in composite parts manufactured by the VARTM process. 
A set of experiments were conducted using a 22-ply, S2 glass, plain weave fabric. All of the 
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preforms were 50.8 cm x 50.8 cm and the distribution media were 45.7 cm x 45.7 cm. The 
panels were manufactured using a line feed source at one end and a line vacuum source at 
the other end. Final part thicknesses were measured at different locations after the parts 
were fully cured. It was found that the final cured thickness of a part changes along the resin 
progression direction. As the distance from the vacuum source increases, the part thickness 
increases, and correspondingly, the fiber volume fraction of the laminates decreases. The 
thickness variation across the panel was close to 6 percent. 

Williams et al. [30] conducted a preliminary experimental study of fabric compression 
during the VARTM process. Nine layers of plain weave E-glass fabric of 660 mm x 360 
mm were injected with unsaturated polyester resin. Two LVDTs, LVDT1 and LVDT2, were 
installed to record the thickness change at the locations of 30 mm and 180 mm away from the 
resin inlet, respectively. Two pressure transducers were positioned at the same locations. 
Four traces were obtained from four nominally identical runs. It was observed that the 
changes in thickness are consistent for the four experiments. Both LVDTs recorded an 
initial, almost instantaneous compaction of the preform as the vacuum was applied, and the 
reduction in thickness continued at a slower rate after the initial change. At LVDT1 location, 
two of the experiments recorded a rapid decrease in thickness at the resin front. This might 
be explained by the initial lubricating effect of the fluid. However, the lubricating effect was 
not apparent in the other two experiments. A possible explanation for the absence of the 
lubrication effect would be that the fabric layers were more closely nested at the beginning 
of the experiment, thus limiting the ability of the fiber network to rearrange. After the 
flow front passed, a steady increase in thickness was observed at LVDT1 location. However, 
when the resin reached LVDT2 location, the decrease in thickness was evident and was 
maintained as the flow continued until the mold filled. After the mold is filled, the resin 
inlet was closed while the vacuum was left on. This resulted in a steady decrease in resin 
pressure, accompanied by a further compaction of the laminate at LVDT1 location, but the 
thickness did not change at LVDT2 location. The mechanisms governing these experimental 
observations are the subject of the further investigation. 
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Detailed knowledge about the response of reinforcing fabrics subjected to compressive 
forces is very important in the investigation of the compaction behavior of the preform dur- 
ing the VARTM process. Kim et al. [31] presented the results of compression and relaxation 
tests using different fiber samples, including E-glass, graphite cloth, mat and unidirectional 
material, and combinations of two different fiber orientations. In the compression test, a nor- 
mal force up to 8.6 MPa was applied to compress the fiber samples at a cross-head speed of 
2 mm/min. The curve of compression force versus sample thickness was measured. Accord- 
ingly, the relation between the stress and the fiber volume fraction was determined. For the 
relaxation test, the machine was stopped after an initial stress was applied to the samples. 
Stress was then measured as a function of time while sample thickness was held constant. It 
was concluded that the compressive response of the reinforcing materials strongly depends on 
the fiber orientation and stacking sequence. When loading and unloading the fiber samples, 
hysteresis effects were observed. Hysteresis may be due to structure relaxation from fiber 
rearrangement or to fiber breakage. Repeating compression on the same sample resulted 
in further compaction of the reinforcement material. It was also shown that the lubricated 
fiber samples were compacted more than the dry reinforcement under the same compres- 
sion force. Effect of the fiber orientation on the stress relaxation was studied. Relaxation 
behavior decreases with fiber alignment. Random mats and 0/90 cloth showed much more 
relaxation than unidirectional fibers. Data of relaxation was fitted with a Maxwell-Wiecher 
viscoelastic model. 

Robitaille and Gauvin [32] gathered the published experimental data for the compaction 
and relaxation of random mats and woven fabrics to identify general trends related to the 
compaction of textile reinforcements used in composite processing. Analytical compaction 
models in the literature were also reviewed, but no published compaction model can fully 
describe the observed behavior of the reinforcements due to the complexity of the problem. 
An alternative to analytical models is to interpolate the experimental data directly. The 
following two equations were used to fit the experimental results obtained from compaction 
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and relaxation experiments: 


( 2 . 1 ) 


V f = AP B 

— = 1 - Ct (1/D) 

Po 

Fiber volume fraction Vf is a function of the compaction pressure P for compaction curves, 
and the ratio of the current compaction pressure P to the initial compaction pressure P 0 
is a function of the time t. for relaxation curves. Eight observed parameters were defined 
based on curve fitting of the experimental data: the initial fiber volume fraction at which 
the compaction pressure P begins to rise (A), the stiffening index of the reinforcements ( B ), 
the average rigidity of the reinforcements (M), the maximum fiber volume fraction reached 
during part production ( Rvj ), the pressure decay after 1 second ( C ) and after 300 seconds 
(P 300 ), and the rate of reduction of the pressure during relaxation ( D ). Thus, comparisons 
of the experimental curves were conducted and the effects of various processing parameters 
were identified. In compaction, as the number of layers rises, the representative rigidity M 
is increased, the initial fiber volume fraction A is increased, and the stiffening index B is 
decreased both for random mats and woven rovings. In addition, the representative fiber 
volume fraction Rvj decreases with a rise of the number of layers for woven materials. For 
successively applied compaction cycles, A is increased, B is decreased, and Rvj is increased 
for both random mats and woven rovings. The representative rigidity M is also increased 
for woven rovings. Compaction curves progressively move to higher V} values. In relaxation, 
the parameter P 30 o is increased (less relaxation) as the applied pressure rises, is decreased as 
the compaction speed rises, and is decreased as the number of layers rises for both random 
mats and woven materials. 

Reorganization of the fiber networks was shown to have a major impact on the com- 
paction behavior of textile reinforcements subjected to successive cycles [33]. Robitaille and 
Gauvin [34] demonstrated the reorganization of the fiber network and the effect of friction. 
The results of sequences of successive compaction cycles applied to dry textiles and to tex- 
tiles saturated in distilled P 2 G and silicon oil were presented. It was shown that the major 
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reorganization mechanism during the first compaction cycles was lateral spreading of the 
rovings. It resulted in a significant rise of maximum fiber volume fraction. The fiber vol- 
ume fraction continued to rise afterwards, but the compaction behavior became essentially 
dependent on the frictional characteristics of the fiber contacts. It was concluded that the 
application of a limited number of compaction cycles prior to the resin injection could be 
used to reduce the volume of the resin rich areas. 

More recently, Chen et al. [35] developed a nonlinear compaction model for fibrous pre- 
forms. A key physical constant, preform bulk compressibility (C&) was defined to establish 
the constitutive law between the applied pressure and the preform bulk volume: 

AV 

ir = ~ C > AP ( 2 - 3 ) 

to 

where, Vo is the initial bulk volume and AV the bulk volume change. A typical pressure- 
thickness curve from the compression test can be divided into three regimes: an initial linear 
regime, a final linear regime, and a nonlinear regime in between. Regime 1 was initiated 
when a minimum amount of pressure was applied to slightly compact the woven fabric lay- 
ers. The fiber cross-sections were brought into closer contact, the bulk deformation of the 
preform came almost entirely from the compressibility of the interstitial space. In the next 
stage, Regime 2, the larger voids had been filled, and the deformation was contributed by 
the compression of both the solid and the deforming voids. It was characterized by the 
nonlinear stage as observed in the experiments. In Regime 3, the cross-sectional geometry 
of the preform has achieved a self-similarity such that the porosity approached a constant. 
A micromechanical conceptual model was proposed to determine the preform bulk com- 
pressibility in these three regimes. The bulk compressibility Cj, is expressed as a function 
of porosity and five parameters: the initial fiber volume fraction, the maximum attainable 
fiber volume fraction, the initial pore volume compressibility, the fiber compressibility, and 
an experimental index. The model prediction was verified by compaction experiments on 
plain-woven fabric preforms and unidirectional non-woven materials. Excellent agreement 

19 


39 



between the theory and experiments had been obtained. 

The research work in reinforcement compaction was mainly carried out for the resin 
transfer molding (RTM) process where the compaction pressure goes as high as 1 MPa, and 
higher compaction rates are involved. Only a few investigations [29,36,37] on the compression 
behavior of fibrous reinforcement during VARTM, i.e., under low compaction pressure (< 1 
atm) were published. 

Hammami and Gebart [36] conducted compaction tests with the maximum compaction 
pressure not exceeding 1 bar. It was found that for low compaction pressures (< 1 bar), 
the number of layers does not have significant influence on the compaction behavior of the 
preform. Two stacking sequences, i.e., [(0°/90°)5] and [(0°)5/(90°)5] were tested to investigate 
the effect of stacking sequence on the compaction response of the reinforcement. The results 
were that the stacking sequence did not make significant difference for the compaction of dry 
reinforcements. However, alternating the reinforcement layers helped achieve a higher fiber 
volume fraction for the compaction of fiber samples saturated with resin. The compressive 
response of the preform was modeled by interpolating the experimental data using dual 
kriging method. A 1-D model was presented to simulate the VARTM process. In this model, 
Darcy’s law and continuity equation were used to model the resin flow, and an equation of 
transverse equilibrium was introduced to take into account the compaction of the preform 
under the vacuum force. Compared with the infusion experiments, the model obtained a 
good estimate of the total filling time for VARTM procedure, but the authors did not present 
any results about the thickness change of the preform. 

Rigas et al. [29] performed several experiments on the compressibility of fiberglass pre- 
forms under low compaction pressures. It was found that by mechanically loading and 
unloading the preform, a maximum fiber volume fraction may be attained prior to applying 
vacuum for resin infusion. The lubrication effect of the resin was also illustrated. In the ex- 
periment, a 16% increase in the fiber volume fraction was observed for the preform saturated 
with resin, compared to the dry fabric under the same load. 
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2.3 Summary 


In our research group, a three-dimensional model, 3DINFIL [38, 39] was developed for the 
simulation of the resin transfer molding (RTM) and resin film infusion (RFI) procedures. 
The physical processes modeled include the flow of a resin through a three-dimensional 
anisotropic preform, heat-transfer in the preform and tooling components, and cure kinetics 
of the resin. The liquid resin was assumed to be an incompressible, Newtonian fluid. The 
governing differential equation for the flow model was based on Darcy’s law and solved by 
the Finite Element/Control Volume (FE/CV) technique. 

Sayre [40] used the computer model 3DINFIL to simulate the isothermal resin flow in 
VARTM process. Experiments were conducted to verify the accuracy of the predicted flow 
front patterns and infiltration times. The flow front patterns were determined to be quali- 
tatively accurate in both the distribution medium and the preform. Simulation results for 
infiltration times over-predicted the experimental measurement by 8 to 10 %. The difference 
may be attributed to an inaccurate porosity and permeability model of the preform used 
in the simulation. In 3DINFIL, the flow model was developed for RTM/RFI process and 
the porosity and permeability were assumed to be constant during the injection procedure. 
However, in VARTM process, the flexible nature of the vacuum bag coupled with the vary- 
ing pressure inside the mold cavity results in a variation of the preform thickness during the 
impregnation. This complex compaction behavior of the preform leads to the variation of 
the porosity and permeability of the preform during infusion process. 

Therefore, the objective of this research was to modify the existing RTM/RFI model to 
include VARTM simulation capabilities and to verify this model with VARTM experiments. 
The new model, 31)1 X I'll. 5.0, features a numerical compaction model incorporated with the 
resin flow model to predict the final thickness and fiber volume fraction of the composite 
part, details of resin flow in the porous preform, and the total infiltration time. A variety 
of element types, including ID rod element, 2D triangle element, 2D quadrilateral element, 
3D tetrahedral element, 3D pentahedral element, and 3D hexahedral element, have been 
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incorporated into the computer code, enabling efficient modeling of the VARTM process in a 
variety of conditions, such as shell geometries, parts with inserts, panels of complex shapes, 
and sandwich structures. 
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Chapter 3 


Model Development 


The VARTM process consists of two important mechanisms: flow of the resin through the 
preform, and compaction and relaxation of the preform under the vacuum force. In addition, 
the resin viscosity and cure kinetics are important for the accurate simulation of the infusion 
process and to determine whether the mold can be filled before the resin gellation. Corre- 
spondingly, flow model, compaction model, and resin model are developed in this chapter to 
simulate the VARTM fabrication procedure of the composites. 


3.1 Flow Model 

The flow model is developed to track the flow of the resin through the distribution medium 
and the preform. Both the high-permeable distribution medium and the preform can be 
modeled as heterogeneous and anisotropic porous media. The resin fluid is assumed to 
be Newtonian (its viscosity is independent of shear rate), and incompressible. Assuming 
that the flow is quasi-steady state, the governing equations for the flow problem are the 
continuity equation for an incompressible fluid, and Darcy’s law of flow through a porous 
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medium [26-29]: 


V • v = 0 


(3.1) 



(3.2) 


where v is the interstitial velocity vector of the resin, q the superficial velocity vector, <p the 
porosity of the porous preform, ^ the viscosity of the fluid, S the permeability tensor of the 
preform, and P r is the resin pressure. Substituting eqn. (3.2) into (3.1) gives the governing 
differential equation of the flow: 

V-(J-VP,.) = 0 (3.3) 

The boundary conditions necessary to solve eqn. (3.3) are: 


• A flow front pressure condition: P rt , , =0 

r ' flow front 

• A constant pressure condition at the inlet: P rinlei = 1 atm 

• The velocity normal to the mold wall is zero: v ■ n = 0, where n is the vector normal 
to the boundary wall 


Note that this is a moving boundary problem. The finite element/control volume (FE/CV) 
method is utilized to track the progression of the flow front. At each time step, the Galerkin 
finite element method is used to solve for the pressure distribution in the fluid. The resin 
velocities are then calculated using eqn. (3.2). With resin velocities obtained, the flow front 
location at that time is determined by means of the control volume technique. The details 
of the FE/CV technique will be addressed in Chapter 4. 
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3.2 Compaction Model 


3.2.1 Compaction Model Formulation 

In addition to the governing equations deriving the flow model, i.e. , Darcy’s law and the 
continuity equation, a new equation is introduced to account for the transverse equilibrium 
inside the mold cavity during impregnation [36]: 

Pc = Pr + Pn (3-4) 

where P c is the total compaction pressure, P r is the resin pressure, and P n is the net pressure 
applied to the preform. For the VARTM process, the external pressure applied is the atmo- 
spheric pressure. It is well accepted that during the flow of the resin in the fiber preforms, 
the total vacuum pressure is shared by the resin pressure and the pressure supported by the 
fiber network [41,42], 

At each time step, once the resin pressure distribution is obtained from the resin flow 
model, the pressure supported by the preform is computed using eqn. 3.4. The normal strain 
in the preform along the transverse direction (e) is a function of the net pressure applied to 
the preform. There is no general constitutive model which can fully describe the response 
of the reinforcement under the compression force. Compaction tests are often conducted to 
determine curves of compressive strain versus compression pressure. With the initial fiber 
volume fraction of the preform (V} 0 ) and thickness of the panel (to) given, the fiber volume 
fraction Vj . displacement ilk and thickness t. could be found by the following equations: 

pi o 

W — £(ly (3.5) 

Jy = o 

t= f (1 + s)dy (3.6) 

Jy = o 
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(3.7) 


V, = V/oj 

The flexible nature of the vacuum bag coupled with the varying pressure inside the mold 
cavity results in the fiber volume fraction and thickness change of the preform during the 
impregnation process. 


3.2.2 Compressibility of the Reinforcement 


Since no general constitutive model is available, the relationship between the compressive 
strain in the preform and the applied pressure is obtained by fitting the compaction test re- 
sults to an empirical model. Two important phenomena are observed during the compaction 
experiments. First, because of the resin lubrication effect, the wet fiber sample saturated 
with resin is compacted more than the dry reinforcement under the same pressure. Second, 
the compressive response of the preform under the compaction force is not elastic, and the 
hysteresis occurs during unloading process [29,43]. 

During the VARTM infusion process, before the resin front approaches, the dry rein- 
forcement is under vacuum compression. Thus, the compressive strain of the preform can 
be calculated from the compaction response of the dry preform during the loading process. 
After the resin passes, the local net pressure applied to the preform decreases as a result of 
the increasing resin pressure. This is equivalent to an unloading process. Accordingly, the 
strain in the wet preform is determined by the compaction response of the resin saturated 
preform during the unloading process. Therefore, the compressive strain varies with the net 
pressure applied to the preform in the form of eqn. 3.8. 


£ = 


fi(P n ) dry compaction , loading 
/ 2 (P n ) wet compaction, unloading 


(3.8) 


For a standard VARTM process, the resin pressure is maintained at 1 atm in the injection 
reservoir and 0 atm at the flow front (vacuum side). Fig. 3.1 presents the two different effects 
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of the flowing resin on the compaction behavior of the preform. In dry condition (3.1a), the 
preform essentially supports the external pressure, P n = P a tm ? and a maximum debulking 
deformation of the dry preform is reached. During infiltration, two deformation mechanisms 
are present in the wet area of the preform: the wetting compaction (3.1b) and the springback 
(3.1c). The wetting compaction is caused by a rearrangement of the fiber network, which is 
created by the lubrication effect of the wetting fluid [30]. Under a given external pressure, 
the lubrication of the dry preform will cause an increase in the preform compaction (e) by an 
additional amount of wetting deformation ( e w ). On the other hand, the presence of the resin 
leads to an increase of the local resin pressure (P r ). According to eqn. 3.4, the net pressure 
applied to the preform P n has to decrease since the total external pressure is constant ( P a tm )• 
Consequently, the preform compaction (e) decreases by an amount of springback deformation 
(ea).The two competing mechanisms, wetting compaction and springback, work together to 
determine the local preform compaction behavior. It is clear that at any time during the 
infiltration process, the net compaction of the preform will depend on the relative magnitude 
of the wetting and springback deformation mechanisms. 


3.3 Resin Model 

The VARTM process uses thermosetting polymeric resins. As the process progresses, the 
resin begins to cure and change viscosity. The VARTM fabrication method is often used to 
manufacture large structures that require long processing times, and the completion of the 
impregnation process before resin gellation becomes a major challenge. Therefore, a model 
is necessary to predict the cure and viscosity change of the resin. 
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Figure 3.1: Compaction mechanisms during VARTM, a) dry compaction, b) wetting com- 
paction and c) springback 
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3.3.1 Cure Kinetics Sub-Model 


A typical expression for the cure kinetics model is given as follows [44]: 

C f t =f(T, a)(l-a)” (3.9) 

where da/clt is defined as the reaction or cure rate, a the degree of cure, T the resin tem- 
perature, /(T, a) a function that depends on the reaction type, and n is the experimentally 
determined constant. The cure rate is a function of resin temperature and degree of cure. 
The VARTM process is often carried out at room temperature. Neglecting the heat liberated 
from the curing resin, the resin temperature is assumed to be constant and equal to the room 
temperature. Differential Scanning Calorimetry and Differential Thermal Analyzer are the 
two common approaches to develop a kinetic model. 

If diffusion of chemical species and convection of the fluid is neglected, the degree of cure 
at each point inside the material can be determined by integrating the expression for the 
cure rate with respect to time in the following manner: 

“ = l O* (3.10) 


3.3.2 Viscosity Sub-Model 

To accurately predict the resin infiltration process, the viscosity of the resin must be known 
as a function of both position and time. Resin viscosity is a complex function of shear rate, 
degree of cure and temperature. There is no existing analytical model which can adequately 
describe this relation. A reasonable approach is to assume that the resin is a Newtonian 
fluid and to measure the viscosity at low shear rates. The measured viscosity can then be 
fit to a mathematical expression relating temperature and time to viscosity. The viscosity 


29 


49 



model used in this study was developed by Castro and Macosko [45]: 


fJ>(T,a) = fi 0 {T) 



-| A(T)+B(T)a 


(3-11) 


where // is the viscosity, /i 0 the viscosity at zero shear rate, T the temperature, a g the degree 
of cure at gel, a the degree of cure, and A and B are parameters which depend on the resin 
temperature. 


3.4 Summary 

In this study, a comprehensive model was developed to simulate the VARTM process. The 
software package, 3DINFIL 5.0, consists of three models: flow model, compaction model, and 
resin model. Fig. 3.2 illustrates the inputs and outputs for each model, and the interaction 
between the three models. 

At each time step, the resin pressure distribution is computed by the flow model and used 
as an input to the compaction model. On the other hand, the output of the compaction 
model, the porosity and permeability of the preform, are updated and influence the resin 
flow in the preform. The viscosity and degree of cure are also updated in the resin model 
at each iteration and passed to the flow model. The viscosity is the required input for the 
resin flow simulation. At each point, the degree of cure (a) is compared with the degree of 
cure at gel (a g ). Once the resin gels (a > a g ), the infusion process can not continue any 
more and the simulation must be stopped. If the resin gels before the mold is filled, the 
process parameters must be adjusted to ensure the complete filling of the mold. The flow 
model, compaction model, and resin model are coupled together as a comprehensive VARTM 
simulation tool to monitor the resin flow pattern in the distribution medium and preform, 
track the thickness change of the preform, predict the infiltration time and final thickness of 
the composite part, and determine whether complete filling can be achieved before the resin 
gels. 
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INPUT MODEL OUTPUT 



Figure 3.2: Inputs and outputs of flow model, compaction model, and resin model 
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Chapter 4 


Finite Element/Control Volume 
Method 


The finite element/control volume (FE/CV) method was selected to solve the resin flow 
problem because of the advantage that this fixed mesh method eliminates the need for 
remeshing the resin-filled domain for each time step, thus the flow simulation for a complex 
geometry can be done rapidly and efficiently. 

In the FE/CV approach, the mold is first discretized into finite elements. By subdividing 
the elements into smaller sub-volumes, a control volume is constructed around each node. 
The concept of fill factor is introduced to monitor the fluid volume in each control volume. 
It is defined as a ratio of the volume of the fluid in the control volume to the total volume of 
the control volume. The fill factor takes on values from 0 to 1, where 0 represents that the 
control volume is empty and 1 when it is filled. The control volumes can be empty, partially 
filled, or completely filled. The numerical flow front is constructed of the nodes that have 
partially filled control volumes. At each time step, fill factors are calculated based on the 
resin velocity and flow into each nodal control volume. If the resin does not gel, the flow 
front is thus updated until all the control volumes are full and the infiltration is complete. 
A flow chart illustrating the numerical procedure is shown in Fig. 4.1. 
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Figure 4.1: Flow chart of the numerical procedure of FE/CV technique 
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4.1 Galerkin Finite Element Formulation 


The governing equation and the boundary conditions for the flow problem are summarized 
as follows. 


V • ( — VP,.) = 0 in O 
P r = P r o n Ti 
v ■ n = 0 on T 2 


where 


Ti UT 2 = 90 and Ti n T 2 = <f> 


(4.1) 

(4.2) 

(4.3) 


(4.4) 


The boundary of 0 is denoted as 90, and Fi and T 2 are complementary parts of 90. 

Following the standard procedure of Galerkin’s method [46], the finite element formula- 
tion of the flow problem can be obtained: 


KijPj = 0 (4.5) 

where Pj is the pressure at each node. K tl is the global resin flow matrix and it can be 
assembled from the element resin flow matrices: 


\Ki,\ = £ [Kj (4.6) 

where 

K. = j ^N m , a N nS dn e (4.7) 

Here, N m and N n are shape functions defined over the element domain O e . 
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The variation of pressure in a finite element is assumed as 


P{x a ) = N m (x a )P m (4.8) 

where, N m is the shape function, and P m is the nodal pressure. After the pressures at all 
the nodes have been solved, the velocities at the centroid of each element is calculated using 
Darcy’s law: 

p (49) 

u a iV m,p- £ m \ ±,u ) 

After the velocity of the resin is obtained, the control volume technique is used to track the 
flow front position. 


4.2 Control Volume Technique 

The control volume technique is a numerical approach for tracking flow font location in fluid 
dynamics problems involving flow with a free-surface [4], Compared to the moving mesh 
methods that track the flow using remeshing schemes, the control volume technique allows 
for approximate tracking of the flow front using a fixed mesh. The first step to implement 
the control volume technique is to form a network of the control volumes, which are formed 
from the finite element mesh. Collectively, the control volumes have the properties of not 
overlapping, and completely filling the mold cavity. The flow front is advanced in this 
technique by calculating the macroscopic fluxes passing from full to empty control volumes. 
The principal advantages of this approach are computational efficiency and stability. 

There are three basic components to the control volume technique: (1) control volume 
formation, (2) flux calculation, and (3) flow front advancement. These are discussed in detail 
below. 
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Figure 4.2: Intra-element notation used in the CV technique for linear quadrilateral elements 

4.2.1 Control Volume Formation 

Each node in the finite element mesh is surrounded by a ‘control volume’ that is composed 
of a collection of sub- volumes from surrounding elements. The size and shape of each control 
volume are determined by the structure of the mesh. 

The mold cavity is first discretized using finite elements, and then each individual element 
is further divided into sub-volumes. Each sub-volume is associated with one of the nodes on 
the element. The number of sub-volumes is equal to the number of nodes in the element. The 
control volume for a particular node is composed of all of the sub-volumes associated with 
that node. An example of the control volume formation in case of quadrilateral elements is 
illustrated in Fig. 4.2 and 4.3. In Fig. 4.2, a 4-noded quadrilateral element is divided into 
four sub-volumes. Points 1, 2, 3 and 4 are the nodes of the element. Point o is the element 
centroid and the midpoints of the four sides of the quadrilateral element are denoted as point 
a, b, c, and d, respectively. The element is divided into four sub-volumes defined respectively 
by the regions, laod , 2boa , 3cob and J^doc. Fig. 4.3 illustrates the formation of the control 
volume for node i by assembling all the sub-volumes associated with the node. 
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Sub-Volume 



Control 

Volume 


Figure 4.3: Assemblage of a control volume for node i from the sub-volumes 

4.2.2 Flux Calculation 

In a manner analogous to the control volume formation, there are two parts to the flux 
calculation. First, the local fluxes into each of the sub-volumes are calculated, then the 
total fluxes into each control volume are calculated by assembling the individual local flux 
contributions. 

At each time step, the finite element method is first utilized to calculate the resin pressure 
and velocity as described in Section 4.1. It is assumed that the velocity of the resin is 
constant throughout each linear element and can be used to compute the local fluxes into 
each sub-volume. Fig. 4.4 illustrates the calculation of the resin flow into a sub-volume for 
a quadrilateral element. Points 1, 2, 3 and 4 are the nodes of the quadrilateral element 
denoted as element ‘e’. The sub- volume associated with node 2 is defined by the region 
a2bo (shaded area). The resin flows into sub- volume a2bo through the internal boundaries 
oa and ob. With the resin velocity v e given, the local flux into the sub- volume, Q e 2 , can be 
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Figure 4.4: Calculation of the local flux into the sub-volume 


calculated as 

Qe 2 = (-v e • ni\ob\t) + ( I'e ■ n 2 \oa\t) (4.10) 

where fq and n 2 are the unit outer normal vectors of the surface ob and oa , respectively, 
and t is the thickness of the element. 

In eqn. 4.10, the first term on the right hand side is the volumetric flow rate across surface 
ob , and the second term is the volumetric flow rate across surface oa. Eqn. 4.10 can be also 
written as 

Qe 2 = -v e ■ {ni\ob\ + n 2 \oa\)t (4.11) 

With a new concept of internal face area vector defined as [11] 

a e2 = {ni\ob\ + n 2 \oa\)t (4.12) 

eqn. 4.11 can be further simplified as 

Qe 2 = -V e ■ a e 2 (4.13) 

Generally, the local flux into each sub-volume associated with node ‘n’ in element ‘e’ can 
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be calculated by 


(4.14) 


U!en ^ en 

Once the local flux into sub-volumes are calculated, the fluxes into each control volume 
can be computed by assembling the individual element contributions: 

m 

= (4.15) 

e=l 

where Q n is the total flux into the control volume associated with node l n\ and m is the 
total number of sub-volumes associated with the node. 

4.2.3 Flow Front Advancement 

A concept of fill factor is introduced to monitor the progression of the flow front [4], The 
fill factor is defined as a ratio of the fluid volume in the control volume to the total volume 
of the control volume. After the flow rates into each control volume have been calculated, 
the fill factors can be updated. Given the current time step (A t), the fill factors from the 
previous step (/*), the calculated flow rate (Q n ), and the volume of each control volume 
( V n ) , the new fill factors can be calculated with 

+ 1) = fl + 4^ (4.16) 

v n 

The fill factor takes on values from 0 to 1, where 0 represents totally empty and 1 totally 
full. When the control volume is partially full, the fill factor is between 0 and 1. The 
numerical flow front is constructed of the nodes that have partially full control volumes as 
shown in Fig. 4.5. The accurate location of the fluid in the control volume is not known, so 
the exact shape of the flow front can not be obtained. The resolution of the flow front is 
determined by the mesh density. 
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Figure 4.5: Actual and numerical flow front 





The time step for the next iteration must be calculated before the solution can proceed. 
The optimal time step would be where the fluid just fills one control volume. If a larger step 
were chosen, the flow front would over-run the control volume and result in a loss of mass 
from the system. The time to fill the partially full control volume ‘n’ is calculated with the 
following relation: 

A t n = (1 ~ /w)ln (4.17) 

Vn 

Once A t n has been calculated for all the partially filled control volumes, the smallest A t n is 
chosen as the time step for the next iteration. 


4.3 Implementation of the Control Volume Technique 

The most important parts to the implementation of the control volume method is to divide 
each individual element into sub-volumes, calculate the volume of each sub-volume, and 
calculate the internal face area vector for each sub-volume. In this section, the details of 
these techniques for the various element types are discussed. 

4.3.1 One-Dimensional Rod Elements 

The notation for rod elements is shown in Fig. 4.6. Points 1 and 2 are the nodes of the rod 
element and o is the element centroid. The two sub-volumes associated with node 1 and 
2 are defined by the line segments lo and 2o. The volume of the two sub- volumes can be 
calculated as 


Eel 

= \lo\S 

(4.18) 

E e2 

= \2o\S 

(4.19) 
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Figure 4.6: Intra-element notation used in the CV technique for linear ID rod elements 



Figure 4.7: Intra-element notation used in the CV technique for linear 2D triangular elements 

where, S is the cross-sectional area of the rod. The internal face area vectors of the two 
sub-volumes are 


CL e l — 

10 5 

IM 

(4.20) 

a e 2 = 

|2o| 6 

(4.21) 


4.3.2 Two-Dimensional Triangular Elements 

The point definitions for triangles are shown in Fig. 4.7. Points 1, 2, and 3, are the nodes 
of the triangle element. Points a, 6, and c are the midpoints of the sides, and point o is the 
element centroid. The three sub-volumes associated with node 1, 2, and 3, are defined by 
the regions, laoc , 2boa , and Scob , respectively. 
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The region laoc can be divided into two triangles, lao and loc by the segment ol. The 
volume of the region laoc is the sum of the volumes of the two triangle regions: 

V e \ = -(|ol x oa\ + |ol x oc\)t (4.22) 

2 

where, t. is the thickness of the element. Similarly, the volumes of the other two sub-volumes 
can be calculated as: 

V e 2 = — ( |o2 x oa | + |o2 x ob\)t (4.23) 

Z 

V e3 = -(|o3 x ob\ + |o3 x oc\)t (4.24) 

As shown in Fig. 4.7, ni, n 2 , and n 3 are three vectors perpendicular to oh , ob. and oc, 
respectively. These vectors have the properties that 



ni\ 

= \oh\ 

(4.25) 


fk\ 

= \ob\ 

(4.26) 


n 3 \ 

= \oc\ 

(4.27) 

Therefore, the internal face area vectors 

; defined for the three sub-volumes are 



= 

(ni - n 3 )t 

(4.28) 

h e 2 

= 

(n 2 - ni)t 

(4.29) 

h e 3 

= 

(n 3 - n 2 )t 

(4.30) 


where, t is the thickness of the element. 
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Figure 4.8: Intra-element notation used in the CV technique for linear 2D quadrilateral 
elements 

4.3.3 Two-Dimensional Quadrilateral Elements 

The point definitions for quadrilaterals are shown in Fig. 4.8. Points 1-4 are the nodes 
of the quadrilateral element. Points a, b, c and d are the midpoints of the sides, and point 
o is the element centroid. The four sub-volumes associated with nodes 1, 2, 3, and 4 are 
defined by the regions: laod , 2boa , Scob , and 4doc, respectively. The volumes of the four 
sub-volumes can be calculated as: 

V e \ = — ( | ol x oh | 4- | oi x od\)t 

Z 

V e 2 = -(|o2 x oa\ + |o2 x ob\)t 

Zj 

I'e 3 = -(|o3 x ob\ + |o3 x oc\)t 
2 

V e 4 = — ( |o4 x oc\ + |o4 x od\)t 

Z 

where, t is the thickness of the element. 

■Hi, n 2 , n 3 , and n 4 are vectors perpendicular to oh, ob, oc, and od, respectively. The 


(4.31) 

(4.32) 

(4.33) 

(4.34) 
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Table 4.1: Definition of sub- volumes for tetrahedral elements 


Sub- Volume Number 

Sub- Volume Vertices 

Internal Surfaces 

1 

lao! c — do 2 00 3 

1,3,4 

2 

2 b 0 i a — eo 4 oo 2 

1,2,5 

3 

3co ! b - fo 3 oo 4 

2,3,6 

4 

4do 3 f - eo 4 oo 2 

4,5,6 


vectors satisfy the following equations: 



N 

= 

o« 

(4.35) 



= 

\ob\ 

(4.36) 


N 

= 

\oc\ 

(4.37) 


N 

= 

\od\ 

(4.38) 

Therefore, the internal face area vectors 

i defined 

for the four sub-volumes are 


CL e \ 

= 

(?1i 

— //| )/ 

(4.39) 

a e 2 

= 

{n 2 

- ni)t 

(4.40) 

a e3 

= 

(m 

- n 2 )t 

(4.41) 

(Le 4 

= 

(n 4 

- n 3 )t 

(4.42) 


where, t is the thickness of the element. 

4.3.4 Three-Dimensional Tetrahedral Elements 

The point definitions for tetrahedrons are shown in Fig. 4.9. Points 1-4 are the nodes of the 
element. Points a - f are the mid-points of the line segments connecting the nodes. Points 
oi _4 are the centroids of the sides of the element, and point o is the element centroid. The 
sub-volumes for the tetrahedron are hexahedrons with vertices defined in Table 4.1. 
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Figure 4.9: Intra-element notation used in the CV technique for linear 3D tetrahedral ele- 
ments 

The six interior integration surfaces s* for the local flux calculation are defined in Ta- 
ble 4.2. The normal vectors of the six surfaces are defined as Hi_g: 


Hi = 

-(002 x oa + oh x 001) 

( 4 . 43 ) 

n 2 = 

-(oo 4 x ob + ob x 001) 
2 

( 4 . 44 ) 

h 3 = 

-(001 x oc, + oc x 003) 

( 4 . 45 ) 

H4 = 

-(003 x ocl + ocl X OO2) 

Z, 

( 4 . 46 ) 

n 5 = 

~{oo 2 x oe + oe x oo 4 ) 

( 4 . 47 ) 

He = 

i(oo 3 x of + of X oo 4 ) 

( 4 . 48 ) 


Note that the magnitude of each normal vector is equal to the area of the surface. The 
internal integral faces associated with each sub-volumes are listed in Table 4.1, and the 
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Table 4.2: Internal integration surface definitions for tetrahedral elements 


Surface Number 

Points 

1 

00 2 CIO i 

2 

oo,/ t boi 

3 

00 3 COi 

4 

oo 3 do 2 

5 

00 4 eo 2 

6 

oo 4 fo 3 


internal face area vectors for these four sub-volumes are calculated as 


ft ei = ni + n 3 + n 4 (4.49) 

ft e 2 = n 2 + n 5 - ni (4.50) 

«e 3 = — n 2 — n 3 — n 6 (4.51) 

a e4 = Hq — n 4 — r ?5 (4.52) 


Look at the sub-volume associated with node 1, the hexahedron lao 3 c — do 2 oo 3 . This 
hexahedron can be divided into three polyhedrons 1 — oo 2 ao i — oo 3 co 3 , and i — oo 3 do 2 . 
These three polyhedrons have the same apex, point 1, and their bases are the internal surfaces 
1, 3, and 4, respectively. Therefore, the volume of sub-volume related with node 1 is the 
sum of the volumes of the three polyhedrons: 

V e i = -n 4 • 1 o + -n 3 • 1 o + -n 4 • 1 o 

1 _ _ . 

= -(n 4 + n 3 + n 4 ) ■ lo 

= ^a e i • lo (4.53) 
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Similarly, 


V e2 = 

1^ 

-a e2 ■ 2 o 

(4.54) 

V e3 = 

1^ 

-« C 3 • 3 0 

(4.55) 

V e4 = 

— fl e 4 • 4o 

(4.56) 


4.3.5 Three-Dimensional Pentahedral Elements 

The point definitions for the pentahedral elements are shown in Fig. 4.10. Points 1-6 are 
the nodes of the pentahedral element. Points ci-i , are the mid-points of the line segments 


connecting the nodes. Points o i _ 5 are the centroids of the sides of the element, and point o 
is the element centroid. The sub-volumes for the pentahedral element are the hexahedrons 
with vertices defined in Table 4.3. The nine interior integration surfaces s* for the local flux 
calculation are defined in Table 4.4. The normal vectors of the nine surfaces are defined as: 


Hi 

= ~{oo 3 x oa + oh x ooi) 

Z. 

(4.57) 


n 2 

= -(ooi x ob + ob x oo 4 ) 

Zj 

(4.58) 


n 3 

= — (ooj x oc. + ob x oo 5 ) 

Z 

(4.59) 


n 4 

= -(005 x od + od x 003 ) 

Z 

(4.60) 


n 5 

= -(oo 3 x oe + o"e x oo 4 ) 

Z 

(4.61) 


n 6 

= ^(004 X of + of X 005 ) 

(4.62) 


n 7 

= ^( 002 X 0 ^ + 0 ^X 003 ) 

(4.63) 


n 8 

= -(oo 2 x oh + o/r x oo 4 ) 

Z 

(4.64) 


n 9 

= -(00 2 x oi + oi x 005 ) 

Z 

(4.65) 
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The magnitude of each normal vector is equal to the area of the surface. The internal 
integral faces associated with each sub-volume are listed in Table 4.3, and the internal face 
area vectors for the six sub- volumes are calculated as 


CL e l ~ 

= ni + 

n 3 + 

Hi 

(4.66) 

@ j e2 ~ 

= n 5 - 

ni - 

fl2 

(4.67) 

a eS ~ 

= n 2 - 

n 3 + 

ne 

(4.68) 

d e 4 - 

= n 7 - 

n,i — 

n 9 

(4.69) 

^e5 

= n 8 - 

n b - 

ft 7 

(4.70) 

^e6 - 

= n 9 - 

~ 

n 8 

(4.71) 


The volumes of the six sub-volumes are calculated as 


f >i 3«ei ' 

Ve2 = 3 ^e2 ' 2 O 

V ‘ s = l S ‘ 3 ' fo 

b l ~ a e4 ’ 40 

f e5 — 3«e5 " 5o 

f eft -«e6 ' 6(1 


(4.72) 

(4.73) 

(4.74) 

(4.75) 

(4.76) 

(4.77) 


4.3.6 Three-Dimensional Hexahedral Elements 

The point definitions for hexahedrons are shown in Fig 4.11. Points 1-8 are the nodes of the 
hexahedral elements, points a-l , are the midpoints of the line segments connecting the nodes, 
points Oi_6 are the centroids of the sides of the element, and point o is the element centroid. 
The sub-volumes are the liexahedrals with vertices defined in Table 4.5. For hexahedral 
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Table 4.3: Definition of sub-volumes for pentahedral elements 


Sub- Volume Number 

Sub- Volume Vertices 

Internal Surfaces 

1 

lao 3 c — do 3 oo 5 

1,3,4 

2 

°2bo 1 a — eo^ oo 3 

1,2,5 

3 

3cOi b - fo 5 oo 

2,3,6 

4 

4go 2 i - do 3 oo 3 

4,7,9 

5 

5ho 2 g — eoj oo 3 

r 5,7,8 

6 

6io 2 h — fo 3 ooj 

6,8,9 


Table 4.4: Internal integration surface definitions for pentahedral elements 


Surface Number 

Points 

1 

oo 3 aoi 

2 

ooj bo 1 

3 

00 5 CO} 

4 

oo 5 do 3 

5 

oo 3 eo^ 

6 

oo 4 fo 5 

7 

oo 2 go 3 

8 

oo 2 hoj 

9 

oo 2 io 5 
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Figure 4.10: Intra-element notation used in the CV technique for linear 3D pentahedral 
elements 
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elements, there are 12 internal integration surfaces needed for the local flux calculation. 
These are defined in Table 4.6. The normal vectors of the 12 surfaces are defined as: 


ni = 

~(oo 3 x oa + oh x ooi) 

(4.78) 

n 2 = 

-(oo 6 x ob + ob x ooi) 

Z. 

(4.79) 

n 3 = 

-(odi x oc + oc x ooi) 

(4.80) 

n 4 = 

— (oo 5 x od + od x 00 i) 

(4.81) 

n 5 = 

-(005 x oe 1 or x 003 ) 

(4.82) 

n 6 = 

i(oo 3 x of + of x oo 6 ) 

(4.83) 

n 7 = 

i(oo 6 xoj + ojx oo 4 ) 

(4.84) 

n 8 = 

-(oo 4 x oil + oh x oo 5 ) 

(4.85) 

rig = 

— (002 x oi + OZ X OO3) 
z 

(4.86) 

nio = 

\{ 002 X 03 + 03 X 00 6 ) 

(4.87) 

nn = 

-(oo 2 x ok + ok x odi) 

(4.88) 

nu = 

-(od 2 x ol + ol X OO 5 ) 

Z 

(4.89) 


The magnitude of each normal vector is equal to the area of the surface. The internal 
integral faces associated with each sub-volume are listed in Table 4.5, and the internal face 
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area vectors for the six sub- volumes are calculated as 


Cl e i — // 1 — r/4 df ng 

«e2 = n 2 — n 1 + n 6 
a e 3 = + n 7 - n 2 

Q e 4 = U4 T '^8 — 'U3 

^e5 = ^9 — ^5 — ^12 

«e6 = nio — ne — n§ 
S e 7 = n 11 — n 7 — nio 
d-e8 = ni2 — '^8 — fill 

The volumes of the eight sub- volumes are calculated as 

f el -Ciel ‘ lo 

f'e2 = -(?e2 ' 2 O 

u, = fa«-& 

Ifi ^ fte4 " 

be5 — gfie5 ‘ 5o 
f e6 g a e6 " 6® 

f e7 = g^e7 ' 7o 
f e8 g ft e8 " 


(4.90) 

(4.91) 

(4.92) 

(4.93) 

(4.94) 

(4.95) 

(4.96) 

(4.97) 


(4.98) 

(4.99) 

(4.100) 

(4.101) 

(4.102) 

(4.103) 

(4.104) 

(4.105) 
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Table 4.5: Definition of sub-volumes for hexahedral elements 


Sub- Volume Number 

Sub- Volume Vertices 

Internal Surfaces 

1 

laoi d — eosoos 

1,4,5 

2 

° 2 boi a — foeoos 

1,2,6 

3 

Scoi b — go 4 oog 

2,3,7 

4 

4 d 0 i c — ho 5 oo 4 

3,4,8 

5 

5io 2 l — £03005 

5,9,12 

6 

6 jo 2 i ~ fo 6 oo 3 

6,9,10 

7 

7ko 2 j - go 4 oo 6 

7,10,11 

8 

8 lo 2 k — hosoo 4 

8,11,12 


Table 4.6: Internal integration surface definitions for hexahedral elements 


Surface Number 

Points 

1 

003C101 

2 

OOgbOi 

3 

00 4 CO 3 

4 

00 5 do ! 

5 

oo^eog 

6 

oosfoe 

7 

oo 6 go 4 

8 

00 4 ho 3 

9 

00 2 IO3 

10 

oo 2 jo 6 

11 

oo 2 ko^ 

12 

oo 2 lo 5 
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Figure 4.11: Intra-element notation used in the CV technique for linear 3D hexahedral 
elements 
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Chapter 5 


VARTM of Flat Panels 


Two composite panels were fabricated by the VARTM process using SAERTEX multi- axial 
warp knit carbon fiber fabric and the SI-ZG-5A epoxy resin. Panel 1 contained one stack 
of the carbon fabric, and panel 2 contained four stacks of the fabric. Each stack of the 
SAERTEX fabric consists of seven layers of unidirectional carbon fibers, which are laid up 
in a quasi-isotropic stacking sequence and knitted together [43]. 

VARTM of the preforms was simulated using the simulation model developed in Chap- 
ter 3 (3DINFIL 5.0). To assess the accuracy of the model, the predicted flow patterns and 
the changes in preform thickness during resin infiltration were compared with the values 
measured in the experiments. 


5.1 Material Characterization 

5.1.1 SI-ZG-5A Epoxy Resin 

SI-ZG-5A epoxy resin is a commercially available anhydride-cure epoxy blend developed for 
VARTM process at A.T.A.R.D laboratories. It features low viscosity at room temperature, 
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which is essential for the VARTM process due to the low injection pressure. 

The cure kinetics model for SI-ZG-5A was obtained using isothermal and dynamic differ- 
ential scanning calorimeter (DSC). Isothermal measurements were made between 60° 6' and 
140°C. A modified auto-catalytic equation was chosen to model the curing reaction for the 
resin [50]: 

da K a m ( 1 — a) n 

(If l _|_ ( < Vi («<•.. • V ' ' 

where 

K = Ae { ~ AE/RT) (5.2) 

The auto-catalytic equation was modified to account for a shift from kinetics to diffusion 
control. In the equation, C is the diffusion constant, ctco is the critical degree of cure at 
T = 0 °K, and acr is the constant to account for the increase in critical resin degree of cure 
with temperature. The experimentally determined values for all the constants in eqns. 5.1 
and 5.2 are listed in Table 5.1. 

A Rheometric Ares System-Five Parallel-Plate Rheometer was used to measure the resin 
viscosity. The resin viscosity was found to fit the following equation: 


ft = 


a n ( A+Ba ) 


a g — a 


(5.3) 


where a g is the gel point degree of cure. The constants used in the model are given in 
Table 5.2. 


Combining equations 5.1 and 5.3, the viscosity of the resin can be predicted for given 
temperatures. Fig. 5.1 shows the viscosity as a function of time at 20°C, 25 0 C, and 30°C'. 
At room temperature, the viscosity of SI-ZG-5A changes less than 5% after 100 minutes. 
Because in the present study, the infiltration procedure can be finished less than 5 minutes, 
the viscosity of the resin is assumed to be constant. The viscosity for SI-ZG-5A was measured 
as 0.34Po • s before the injection. 
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Table 5.1: Parameters used in SI-ZG-5A kinetics model 


Parameter 

Value 

Unit 

A E 

61.0 

KJ j mole 

A 


s~ x 

m 

0.551 


n 

1.0 


C 

40 


«co 

0.216 


OiCT 


K~ l 


Table 5.2: Parameters used in SI-ZG-5A viscosity model 


Parameter 

Value 

Unit 

E, 

54.8 

KJ /mole 

\ 

8.3 x 10“ n 

Pa ■ s 

A 

4 


B 

6.25 


OLg 

0.60 
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5.1.2 SAERTEX Multi-Axial Warp Knit Carbon Fiber Preform 


The permeability of the SAERTEX Multi-Axial Warp Knit (MAWK) carbon preform was 
measured in the three mutually perpendicular principal material directions. These directions 
are in-plane along knitting, in-plane normal to knitting, and through the thickness (TTT) 
direction. The power law equation was chosen to fit the permeability data as a function of 
the fiber volume fraction: 

A = a(Vf) b (5.4) 

where, S is the permeability in m 2 , Vf is the fiber volume fraction, and a and b are fit 
constants. The fit constants are presented in Table 5.3, and Fig. 5.2 shows the plots of the 
permeabilities versus fiber volume fraction. 

Compaction tests were conducted to develop an understanding of the compressive re- 
sponse of the SAERTEX fabric under the low compaction pressure during VARTM process- 
ing. Both dry and wet compaction tests were performed. For the wet preform test, the 
preform was initially impregnated with SI-ZG-5A resin. The test results of the compaction 
pressure versus the resulting strain are shown in Figs. 5.3 and 5.4. Note that test data were 
not recorded during the loading portion of the wet compaction experiment for the 4-stack 
fabric preform. Because of the resin lubrication effect, the compression response of the wet 
fibrous preform was found to be quite different from that of the dry fiber mat. The hystere- 
sis phenomenon was observed during the unloading process. During the VARTM infusion 
process, before the flow front passes by, the dry reinforcement is under vacuum compression 
and the strain of the preform can be calculated from the compaction response of the dry 
preform during the loading process. After the resin passes by, the local net pressure applied 
to the preform decreases as a result of increasing resin pressure. This is equivalent to an 
unloading process. Therefore, the unloading portion of the wet compaction test essentially 
simulates the actual VARTM process but in a quasi-static condition. Empirical equations of 
preform strain as a function of pressure for unstitched MAWK fabric are obtained by fitting 
the compaction test data to the mathematical expression shown in eqn. 5.5 for 1-stack of 
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Table 5.3: Fit constants used in unstitched SAERTEX carbon fabric permeability model 



a 

b 

In-plane along knitting 

2.0 x 10- 13 

-7.62 

In-plane normal to knitting 

4.0 x 10- 13 

-5.36 

TTT 

6.0 x 10“ 15 

-7.29 


fabric and eqn. 5.6 for 4-stacks of fabric: 


t = 


0.22(1 — °- 026P «)) dry compaction, loading 

0.028 + 0.20 1 wet. compaction, unloading 


(5.5) 


£ = 


0.24(1 — °- 025P «)) dry compaction, loading 


0.082 + 0.24, 


wet compaction , unloading 


(5.6) 


3.82+P n 

where, e is the strain of the preform, and P n is the net pressure applied to the preform. 


Fig. 5.5 and Fig. 5.6 show the compaction pressure versus the corresponding fiber volume 
fraction for the dry and wet preforms. The curves shown correspond to the loading portion 
of the dry compaction test and unloading portion of the wet compaction test. The maximum 
achievable fiber volume fraction of the dry preform under full vacuum was 54% and 52% for 
1- stack and 4-stack fabric, respectively. However, the wet preforms could be compacted to 
the fiber volume fraction of 58%. This result illustrates the significant lubrication effect of 
the wetting fluid on the compaction behavior of the preform as discussed in Chapter 3. 


5.1.3 Distribution Medium 

The distribution medium used in this study contains three layers of nylon mesh screen. 
The nylon screen is in-plane isotropic. The in-plane and transverse permeabilities of the 
distribution medium were measured to be 8.31 E — 09 m 2 and 5.491? — 10 m 2 , respectively. 
The porosity of the distribution medium is 0.75. 
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Figure 5.5: Fiber volume fraction achieved for 1-stack S VERTEX fabric in the dry and wet 
compaction tests 
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Figure 5.6: Fiber volume fraction achieved for 4-stack SAERTEX fabric in the dry and wet 
compaction tests 

5.2 Experiments 

Two flat panels (60.96 cm by 30.48 cm) were fabricated by the VARTM process using SAER- 
TEX multi-axial warp knit carbon fiber fabric and the SI-ZG-5A epoxy resin [43]. Panel 
1 contained one stack of the fabric while Panel 2 contained four stacks of the fabric. The 
initial preform thickness was 0.19 cm for Panel 1 and 0.82 cm for Panel 2. 

A 1.27 cm thick aluminum tool was instrumented to allow for measurement of pressure 
and displacement of the preform during the VARTM infiltration. The pressure sensors and 
Linear Variable Displacement Transducers (LVDT) were installed to monitor the pressure 
and thickness change of the panels at different locations (see Fig. 5.7). The pressure sensors, 
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labeled PI through P 3 were mounted at the tool surface beneath the fiber preform. PI was 
2.5 cm from the preform edge on the resin inlet side, P 2 was at the center of the preform, 
and P 3 was 2.5 cm from the preform edge on the vacuum side. LVDTs, labeled D 1 through 
D3. were supported above the vacuum bagged preform by a rigid beam. D 1 was 3.9 cm from 
the preform edge on the resin inlet side, D2 was at the center of the preform, and D 3 was 3.9 
cm from the preform edge on the vacuum side. The preform was cut and placed on the tool. 
A layer of Armalon fabric was placed above the preform to serve as the release film. The 
distribution medium was laid on top of the preform and release layer. The medium was cut 
smaller than the preform allowing a 1.27 cm gap between the edge of the medium and the 
edge of the preform along the length. The distribution medium ended at a distance of 2.54 
cm before the end of the preform. These gaps prevent race-tracking of the resin as it flows 
through the media and the preform. A 5.08 cm portion of the medium was set off of the panel 
at the inlet side for inlet tube placement (Fig. 5.7). The outlet tube located on the vacuum 
side of the part was connected to a resin trap and vacuum pump. The bagging procedure 
was completed when the preform, medium and tubing were sealed to the aluminum tool 
using a conformable vacuum bag and sealant tape. 

SI-ZG-5A resin was degassed at room temperature under full vacuum for approximately 
one hour prior to injection. Once the system was equilibrated and all air leaks were elimi- 
nated, the resin was allowed to flow into the distribution medium and fiber preform. The 
flow front at the top surface of the part was recorded using a digital video camcorder. The 
bottom flow front position was obtained from the tool mounted pressure sensor responses 
recorded during the test. Typically, the presence of the flow front at the tool surface is 
indicated by a decrease in the vacuum measured by the sensor. When the preform was fully 
impregnated, as evidenced through visual inspection, the inlet and then the outlet tube were 
clamped. The part was then placed in an oven to cure. 
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Distribution Medium 




Preform 


Figure 5.7: The VARTM instrumented system [43] 
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5.3 Simulation 


In this study, except the side edges of the panels where there is no distribution medium, the 
resin flow condition is uniform across the width of the panels. Therefore, the resin velocity 
in the width direction is negligible and the resin progression during the VARTM process is 
reduced to a two-dimensional flow problem. 

The design of the infusion scheme in this study is worthy of special attention. As shown in 
Fig. 5.8 (a), the distribution medium is set off of the panel at the inlet side and wraps around 
the left edge of the panel. Therefore, the resin can flow into the preform from two directions, 
from above through the top surface and from the left edge of the preform. However, in the 
experiment, because no clamping system was used to ensure contact between the medium 
and preform around the sharp corners of the preform, the distribution medium could be 
detached from the preform at the left edge as shown in Fig. 5.8 (b). In this study, both cases 
(Fig. 5.8 (a) and (b)) were considered and two different models were created and named as 
Model I and Model I, correspondingly. Another model as shown in Fig. 5.8 (c) was also 
created to compare the effect of the different infusion schemes. This model is denoted as 
Model HI. In model HI, no portion of the distribution medium is set off of the panel, and the 
resin tube is placed above the preform at the left end. Interesting simulation results of the 
three models will be discussed in this section. 
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(a) 



(b) 



(c) 

Figure 5.8: Three different VARTM infusion schemes, (a) Model I (b) Model I (c) Model HI 
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5.3.1 Panel 1 


Shown in Figures 5.9, 5.10, and 5.11 are comparisons between the calculated and measured 
infiltration times versus corresponding top and bottom flow front positions for model I, model 
E, and model IE, respectively. The simulation results of model E agree with the experimental 
data very well while model I overpredicted the infiltration time by approximately 1.7 times. 
This indicates that in the experiment, the distribution medium did not perfectly conform to 
the shape of the preform at the left end and model E precisely simulated the actual VARTM 
lay-up. From Fig. 5.10, the flow front on the top surface leads the resin front on the bottom 
surface before the end of the distribution medium. However, after the end of the distribution 
medium, both the top and bottom flow fronts slow significantly, and the bottom flow front 
catches up with the top resin front. This illustrates the importance of the distribution 
medium in reducing the infiltration time. Initially one might think that the infusion scheme 
in model I could shorten the infiltration time because the resin enters the preform through 
the distribution medium from above and the left end simultaneously. However, this results 
in the competition for the resin between the flow along the distribution medium ( Q \ ) and 
the flow into the preform from the left end (Q 2 ) as shown in Fig. 5.12. As shown in Fig. 5.9, 
this competition delays the resin infiltration into the preform by 90 seconds and leads to a 
significant increase of the total infiltration time compared to the infusion scheme in Model E. 
The resin infusion scheme in Model IE was shown to require least time to fill the panel. The 
total filling time predicted by model E and model El is 151 and 140 seconds, respectively. 

Fig. 5.13 shows the displacements measured by the three LVDTs. Fig. 5.14 presents the 
displacements at these three locations calculated by model E, which was shown to accurately 
predict the resin infiltration. The measured and calculated displacements are compared 
in Figures 5.15-5.17. The calculated thickness change of the preform agrees well with the 
experimental measurements qualitatively. Before the resin inlet opens, the dry preform is 
compacted under vacuum and initial displacements are induced. After the infusion process 
begins, the presence of the resin affects the compaction of the preform by two different 
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Figure 5.9: Flow front location comparison between Model I prediction and experiment 
measurement for Panel 1 
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Figure 5.10: Flow front location comparison between Model II prediction and experiment 
measurement for Panel 1 
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Figure 5.11: Flow front location comparison between Model HI prediction and experiment 
measurement for Panel 1 
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Figure 5.12: Competition for resin during VARTM with the infusion scheme in Model I 

mechanisms [43]. First, the lubrication effect of the resin causes rearrangement of the fiber 
network and an increase in the preform compaction. This mechanism is called the wetting 
compaction effect of the resin. On the other hand, the increase of the resin pressure leads to 
a decrease of the pressure applied to the preform. Consequently, the amount of the preform 
compaction decreases. This is called the springback mechanism. During the infiltration 
process, the net compaction of the preform depends on the relative magnitude of the wetting 
and springback deformation mechanisms. Therefore, both the simulation and the experiment 
find that the compaction responses of the preform at the three LVDT positions are different. 
D1 moderately decreased after the resin passed by, while D2 increased after the flow front 
approached. The model also predicted a temporary increase of the displacement D3 due 
to wetting compaction, although it was not observed in the experiment. After the flow 
front reached the end of the distribution medium, Dl, D2, and D3 all decreased rapidly due 
to the sharp increase in the resin pressure. The sudden increase in preform thickness was 
observed at 115 seconds in the experiment, and the phenomenon was predicted to occur at 
116 seconds in the simulation. Both the simulation and the experiment find that at the end 
of the infiltration process, Dl < D2 < D 3. This indicates that the panel is less compacted 
on the resin inlet side and more compacted on the vacuum side. 
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Figure 5.13: Measured displacements of Panel 1 

The computed displacements do not match the measured value quantitatively. There are 
three reasons to explain about the differences. First, the initial thickness of the preform is 
not uniform. This also explains the difference between the measured initial displacements 
of the dry preform under the full vacuum at the three LVDT locations. The value of the 
initial thickness of the preform used in the simulation is the averaged thickness of the panel 
and this induces error. Second, the compressibility of the preform (i.e., the equations of 
the compressive strain as a function of the pressure) used in the simulation was measured 
under quasi-static conditions, which may not truly represent the actual dynamic VARTM 
compaction conditions. Also, the error in fitting the compaction test results contributes to 
the disagreements. 
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Figure 5.14: Predicted displacements of Panel 1 














5.3.2 Panel 2 


As discussed in Section 5.3.1, model E shown in Fig. 5.8 (b) best represents the actual 
VARTM lay-up in the experiment, and is thereby adopted to simulate the infiltration of 
Panel 2, which contains 4 stacks of the SAERTEX carbon fabric. 

The computed and measured infiltration times of the top and bottom resin front positions 
are compared in Fig. 5.18. The flow pattern predicted by the model agrees well with the 
measured values on the top surface before the end of the distribution medium. However, on 
the bottom surface, the predicted resin infiltration was much slower than measured. The 
reason for the disagreement is that for the unstitched preform containing more than one stack 
of the fabric, the permeability model used in the simulation is not accurate enough. The 
current permeability model is obtained by the permeability tests, and the mechanical load 
frame used in the tests can not hold a low compaction load. As a result, the permeability 
of the MAWK fabric can only be measured at a 45% or higher fiber volume fraction. The 
permeability of the preform with a lower fiber content was extrapolated from data points 
for 45% and higher fiber volume fractions. However, for preforms containing more than one 
stack of the fabric, when the pressure applied becomes low and the fiber volume fraction of 
the preform drops below 45%, the gaps between the stacks could increase significantly and 
form flow pathways for the resin at the interply interfaces. That is, the effective permeability 
of the whole structure could increase drastically when the fiber content is below 45%, and the 
empirical permeability model obtained from fiber volume fraction data 45% or higher does 
not apply any more. Note that the fiber volume fraction of the uncompressed 4-stack fabric 
preform is 40%. During the infiltration process, the pressure applied to the preform keeps 
decreasing as a result of increasing local resin pressure. Therefore, the compaction pressure 
becomes very low and the fiber volume fraction of the preform is then likely to be in the 
range between 40% and 45%, where the current permeability model can not predict the 
permeability of the preform accurately, which leads to the inaccurate prediction for the resin 
flow pattern in the preform. The 1 -stack fabric preform does not have this interply interface 
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flow problem. Hence the permeability model works well in the simulation of VARTM process 
for Panel 1. 

The measured and predicted transverse displacements of Panel 2 at the three LVDT 
positions are presented in Fig. 5.19 and Fig. 5.20, respectively. The predicted thickness 
change agrees well with the experiment measurement qualitatively. For Panel 2, at all the 
three locations, the wetting compaction mechanism dominated over the springback mecha- 
nism when the resin front approached and wetted out the preform. However, the local resin 
pressure kept increasing with the progression of the flow front. This led to the persistent 
decrease of the net pressure applied to the preform. After the maximum displacements were 
reached, the displacements decreased as a result of the significantly decreased compression 
force. Both the simulation and the experimental results show that, like Panel 1, Panel 2 is 
also less compacted on the resin inlet side and more compacted on the vacuum side. 
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Figure 5.18: Flow front comparison between simulation and experiment for Panel 2 
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Figure 5.20: Predicted displacements of Panel 2 





Chapter 6 


VARTM of Sandwich Structures 


Sandwich structures are being considered for use in primary aircraft structures due to their 
light weight and good resistance to shearing. NASA Langley Research Center is investigating 
a new form of sandwich structures that contain tli rough-th c-tli i ckness reinforcement. Recent 
studies [47] have shown that transverse stitchings can improve damage tolerance of these 
structures, making them better suited for service in structural applications. It is envisioned 
that the stitched foam core preform could be resin infused by the VARTM process. 

The primary objective in developing the VARTM process of stitched foam core sandwich 
structures is the complete filling of the preform with adequate wetting of the fibrous media. 
Thus, it is very important to investigate the details of the resin flow through the stitched 
foam core preforms. The VARTM simulation model developed in Chapter 3 (3DINFIL 5.0) 
was employed to simulate the flow field during VARTM of the stitched foam core panels. 
The simulation results were verified by experiments. The effect of the stitching density on 
the filling time and infiltration pattern was examined. 
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6.1 Material Characterization 


The stitched foam core composite panels consisted of SAERTEX carbon fiber multiaxial 
warp knit face sheets attached to a Rohacell 31 IG polymethacrvlimide (PMI) foam core. 
The face sheets and foam core were stitched together in the transverse direction using a 1600 
denier Kevlar/PVA thread. The resin used in this study was a room temperature cure Dow 
Derakane 510A-40 vinyl ester. A plain weave nylon screen was used as the high-permeable 
distribution medium. The porosity and permeability of the carbon fabric and nylon screen, 
and the resin viscosity are required inputs for the computer model. 

6.1.1 Porosity and Permeability 

The SAERTEX multi axial warp knit carbon fabric contains seven layers of unidirectional 
carbon fibers laid up in a quasi-isotropic stacking sequence. The seven layers are knitted 
together, and the knitted unit is referred to as a ‘stack’. The permeability of the fabric 
was measured using the steady-state measurement technique in a previous study [48]. The 
fiber samples used in the permeability tests contain 8 stacks of the carbon fabric, which 
were stitched together through the thickness using a modified lock stitch and Kevlar thread. 
The stitch rows were 5.08 mm apart, and the stitch step was 3.75 mm. Assuming that the 
material has orthotropic flow properties, permeability measurements were made in the three 
mutually perpendicular principal directions. For the stitched carbon fabric, these directions 
are in-plane parallel to the stitching, in-plane normal to the stitching, and transverse to the 
stitching (TTT) direction. The permeability of the carbon fabric varies with fiber volume 
fraction. The power law equation was used to fit the permeability data in each direction: 

S = a(V f ) b (6.1) 

where S is the permeability in m 2 , Vf is the fiber volume fraction, and a and b are fit 
constants. The fit constants are listed in Table 6.1, and Fig. 6.1 shows the power-law fit 
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Table 6.1: Fit constants in the permeability model for stitched SAERTEX carbon fabric [48] 



a 

b 

to stitching 

9.03 x 10“ 15 

-12.18 

T to stitching 

3.27 x 10“ 15 

-12.53 

TTT 

1.26 x 10 -14 

-7.95 


curves for the permeability parallel, normal, and transverse to the stitching as a function of 
fiber volume fraction. 

The permeability model was coded into the material database of the computer simulation 
package, 3DINFIL 5.0. With given fiber volume fraction, the permeability of the carbon face 
sheets could be calculated by the code. In this study, the fiber volume fraction (V)) of the 
carbon face sheets were assumed to be constant and took the value of 0.53. Accordingly, the 
porosity ( 0 ) of the carbon fabric was calculated to be 0.47 using the following equation: 

0 = 1- V f (6.2) 

Sayre [40] measured the porosity and permeability of the nylon screen used as the dis- 
tribution medium. By examining the geometry of the distribution medium and calculating 
the ratio of the open cross-sectional area to the total cross-sectional area, the porosity of 
the nylon screen was determined to be 0.51. In the permeability measurement, a radial 
flow, advancing front technique with a center-port injection was first used to determine the 
orientation of the in-plane principal directions for the permeability tensor. The circular flow 
patterns were observed in the radial infiltration test and this indicated that the nylon screen 
is in-plane isotropic. One-dimensional advancing front technique was then used to measure 
the in-plane permeability of the distribution medium, which was found to be 2.92 E — 09 m 2 . 

The nylon screen is very thin and thus, the transverse direction fills immediately as 
the fluid flows in the in-plane directions. This condition prohibits that any realistic trans- 
verse permeability measurements be performed. Therefore, the transverse permeability was 
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Figure 6.1: Permeability of stitched SAERTEX carbon fabric 




arbitrarily assumed to be 1.0 E — 08 rrr [40]. 


6.1.2 Resin Viscosity 

As described in Chapter 3, the resin viscosity (p) is often characterized as a function of 
temperature (T) and degree of cure {a): 

V = f(T,a ) (6.3) 

The degree of cure {a) of the resin varies with temperature (T) and time (t): 

a = g{T,t) (6.4) 

Substituting eqn. 6.4 into eqn. 6.3 yields 

I 1 = f(T,g(T,t)) (6.5) 

Therefore, the resin viscosity (p) can be modeled as a function of temperature (T) and time 

(*): 

A' = p(T, t) (6.6) 

Because injection of the sandwich panels with 510A-40 vinyl ester was performed at room 
temperature, the resin temperature is assumed to be constant during the VARTM process. 
Thus, the resin viscosity model used in the study can be further simplified and takes the 
form of the following relation: 

p,. = p,.(£) (6.7) 

where p r is the resin viscosity at room temperature. 

Experiments were conducted to determine the viscosity model of the Dow Derakane 
■510A-40 vinyl ester. The viscosity measurements were made with a Brookfield Model DV- 
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Figure 6.2: Viscosity trace of Derakane 510A-40 resin with 0.2% CoNap and 1.0% MEKP 

HI Concentric Cylinder Digital Rheometer. Isothermal measurements at room temperature 
were conducted. The vinyl ester resin was mixed with a cobalt napthenate catalyst at 0.2 
weight % and a methyl ethyl ketone peroxide (MEKP) promoter at 1.0 weight %. Three 
different batches of the vinyl ester were used in the VARTM of the sandwich structures. 
Each has a subtle difference in the resin viscosity, and thus the viscosity of each batch was 
evaluated. Fig. 6.2 shows the measured viscosity as a function of time for the three batches 
of 510A-40 vinyl ester resin. The initial viscosities for the three batches were 0.475 Pa-s, 
0.425 Pa-s, and 0.4 Pa-s, respectively. The resin gelled after approximately 45 minutes. To 
avoid incomplete impregnation of the preform, the mold must be filled before the resin gels. 
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6.2 Experiment 


A dry, stitched foam core preform consists of two 3mm thick carbon face sheets attached to a 
polymerthacrylimide (PMI) foam core by stitching through the thickness of the preform with 
a 1600 denier Kevlar/PVA thread. Three preforms with different stitch row spacings were 
resin injected to evaluate the effect of the stitch density on the filling time and infiltration 
pattern. These included stitch row spacings of 6.25 mm, 12.7 mm, and 25.4 mm. The pitch 
was held constant at 4 stitches per 25.4 mm. All the preforms were approximately 30.5 cm 
wide and 64.5 cm long and used a 12.7 mm thick foam core. 

A VARTM process for resin infiltration of stitched foam core preforms was developed. 
The stitched foam core preform was placed on a glass tool plate. The tool plate was treated 
with a release agent over the area on which manufacturing took place. The glass plate was set 
on a transparent polycarbonate table. The glass tool plate and transparent polycarbonate 
table were selected in order to observe resin infiltration into the lower face sheet of the 
preform. A porous peel ply was placed directly on top of the preform. Then, a nylon 
high-permeable distribution medium was placed on top of the peel ply and aligned with the 
injection edge of the panel. The distribution medium was cut to leave approximately a 25.4 
mm space between the edge of the medium and the vacuum edge of the preform and 12.7 
mm from the sides of the medium and the sides of the preform. A resin distribution tube 
was placed on top of the distribution medium at a distance of 25.4 mm from the edge. The 
resin distribution tube was prepared by drilling 3.18 mm diameter holes 25.4 mm apart. The 
number of holes was determined by the width of the panel. One end of the tube was plugged 
with a brass fitting. A typical lay up of the VARTM experiment is shown in Fig. 6.3. 

After the lay-up was completed, glass bleeder/breather material was cut to the same 
width as the preform and placed on the tool plate at the end opposite to the resin tube. 
The entire assembly was then sealed in a vacuum bag. A vacuum port was centered on the 
bleeder/breather materials. The other end of the resin tube was plugged and a vacuum was 
applied to expel any air. The entire assembly is shown in Fig. 6.4. 
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Figure 6.3: Schematic diagram of the VARTM lay-up for the infiltration of the stitched foam 
core sandwich panel 



to resin supply 


Figure 6.4: VARTM processing set-up for a stitched foam core composite panel 
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Derakane 510A-40 vinyl ester resin was used in this study. Approximately 1300g resin was 
used to inject the 30.5 cm x 64.5 cm panels. The vinyl ester resin with 1.0 % MEIvP and 0.2 % 
cobalt napthanate has a processing time of approximately 45 minutes at room temperature. 
Prior to injection, the resin was allowed to degas. Once the system was equilibrated and all 
air leaks were eliminated, the resin tube was placed into the resin source and the infusion 
process began. After the preform was completely infiltrated, the resin distribution tube was 
clamped off while the vacuum source was left on until the part completely gelled. The panel 
was then cured at room temperature for 24 hours before removing from the mold. 

Flow visualization experiments were conducted to record the flow front positions in the 
upper and lower carbon face sheets. In order to measure the flow front position in the upper 
and lower carbon face sheets, straight lines were drawn at measured intervals on the surface 
of the vacuum bag and on the bottom surface of the glass tool plate, respectively. The 
complete preform layup with grid markings is shown in Fig. 6.5. Visual observations were 
made in two ways. On the top surface, a digital camera was used to capture images of the 
flow front as a function of time. A minute/second timer was used to record the infiltration 
time. The timer was started when the resin filled the resin distribution tube. Then, as the 
flow reached a grid mark, the time was recorded and a digital image was captured. Using an 
angled mirror, the images from the bottom surface were reflected and recorded using a video 
camera. Subsequent flow front position versus time readings was determined from review of 
the videotape. 


6.3 Simulation Model 

The resin infiltrates the preform in two ways. First, the resin moves across the liigh- 
permeable distribution medium, infiltrates the top face sheet, flows through the non-porous 
foam by ‘leaking’ through the holes created by the needle penetrations of the stitching, and 
finally infiltrates the bottom face sheet. Therefore, in order to simulate resin infiltration of 
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Figure 6.5: Experimental setup of a stitched preform. The grid marks are used to determine 
the flow front location as a function of time 

stitched foam core sandwich structures, a finite element mesh was created including a high- 
permeable distribution medium, and two porous multiaxial warp knit face sheets connected 
by the transverse stitches. The closed-cell nature of the PMI foam core makes it impervious 
to fluid, thus the foam core is excluded from the resin flow model. 

6.3.1 Model of the Resin Flow along the Transverse Stitching 

Resin leaks through the needle penetrations of the stitching to infiltrate the preform in the 
transverse direction. The resin path along the stitches from the upper face sheet to the 
lower face sheet can be seen in a cross-sectional view of the stitching in Fig. 6.6. Therefore, 
the resin flow along each transverse stitch can be regarded as one-dimensional flow and 1- 
D rod elements can be used to model the stitches. This type of finite element model was 
named as ‘the channel model’, where the needle penetrations in the foam core were described 
by a series of channels that permit 1-D resin flow from the upper to the lower face sheet. 
However, the finite element mesh for the channel model is very large and results in expensive 
computational coast. For example, a panel of the dimension 30.5 cm x 64.5 cm with the 
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stitch row spacing of 6.25 mm and the pitch of 4 stitches per 25.4 mm contains about 4365 
transverse stitches, and requires more than 10000 nodes in the whole finite element mesh. 
Therefore, a ‘strip model’ was proposed to simulate the resin flow in the stitched sandwich 
panel more efficiently. In the strip model, instead of modeling the single stitches using the 
1-D rod elements, a two-dimensional strip encompassing the transverse stitching was used to 
represent one row of stitching. The channel model and strip model are compared in Fig. 6.7. 
By setting the permeability normal to the transverse stitching threads to be zero, these strips 
allow resin flow in the transverse direction of the foam core only. 



Figure 6.6: Cross sectional view of stitching, resin flow from top face sheet to bottom face 
sheet along the stitching 

The strips have the same dimensions of 61.9 cm x 1.27cm. The thickness of the strips was 
assumed to be 1.60 mm, which is the diameter of the holes created by the needle penetrations 
of the stitching. With one strip representing one row of stitching, the panel with the 6.35 
mm stitch spacing has a total of 44 strips, the panel with the 12.7 mm stitch spacing has a 
total of 23 strips, and the panel with the 25.4 mm stitch spacing has a total of 12 strips. 

The porosity of each strip is determined by considering the ratio of the pore area to the 
total area of the strip as follows: 

cf> = (6.8) 

strip 
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Figure 6.7: Schematic of the channel model and strip model 

area of the pores is calculated by considering the area of the hole made by the needle 
subtracting the area of the stitching threads through the holes: 


Apores — X (Ahoie Afh rea( j^ (6-9) 

where, n is the total number of transverse stitches encompassed in one strip, A^oie is the 
area of the needle penetration and A thread is the area of the reinforcing kevlar thread. Using 
equations 6.8 and 6.9, the porosity of the strips was determined to be about 0.18. 

It is difficult to measure the permeability of the porous region (i.e. , strips) in the foam core 
by experiments. Therefore, a theoretical model, straight capillary model [49] was modified 
to estimate the transverse permeability of the strip. The model is illustrated in Fig. 6.8. 
The porous strip is represented by a bundle of straight parallel ‘capillary unit’. The capillary 
unit is a circular hole holding a concentric cylinder in it as shown in Fig. 6.9. The radius 
of the cylinder, r } , is less than the radius of the hole, r 2 . The circular hole represents 
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Figure 6.8: Modified straight capillaric model 


the hole created by the needle penetration of the stitching and the cylinder represents the 
thread through the hole. The flow 0 through a capillary unit is then given by solving the 
Navi ('r- Stokes equations: 


Q 


s i 

//, cl y 


( 6 . 10 ) 


where 


5 


Ci and C 2 are defined as 


p-2 i 

27 r / (-r 3 + C\r ln r + C 2 r)dr 

J n ^ 


( 6 . 11 ) 


<?! 

C 2 


4(lnr 2 — lnri) 
ri 2 In r 2 — /’ 2 2 In r\ 
4(ln r\ — In r 2 ) 


( 6 . 12 ) 

(6.13) 


ji is the viscosity and ^ the pressure gradient along the capillary unit. 

If there are n such capillary units in the strip, the flow per unit area of cross-section of 
the strip, q will be 

nS 1 dP 
Astrip d dy 


(6.14) 
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Figure 6.9: A capillary unit 
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As the flow can also be expressed by Darcy’s law 


SdP 
// dy 


(6.15) 


It follows that 


S 


nS 


A 


strip 


(6.16) 


The permeability of the strip in transverse direction was calculated to be 4.02 E — 9 m 2 . 


6.3.2 Finite Element Mesh 


The finite element meshes for the three sandwich panels are shown in Fig. 6.10, Fig. 6.11, 
and Fig. 6.12, respectively. 3-D hexahedral elements were used to discretize the distribution 
medium and carbon face sheets, and 2-D quadrilateral elements were used for the simulation 
of resin flow in the porous region (i.e. , strips) of the foam core. Table 6.2 lists number of 
elements and nodes used in the simulation model. 


Table 6.2: Number of elements and nodes used in the simulation models 


Stitch row 
Spacing (mm) 

Number of 
2-D elements 

Number of 
3-D elements 

Total number 
of elements 

Number of 
nodes 

6.35 

2160 

3744 

5904 

7934 

12.7 

1656 

2600 

4256 

5558 

25.4 

864 

2032 

2896 

4057 
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Figure 6.10: Finite element mesh for the stitched foam core preform with 6.35 mm stitch 
row spacing 
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Figure 6.11: Finite element mesh for the stitched foam core preform with 12.7 mm stitch 
row spacing- 
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Figure 6.12: Finite element mesh for the stitched foam core preform with 25.4 mm stitch 
row spacing 
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6.4 Results 


Three stitched preforms with different stitch row spacings were resin injected to examine 
the effect of the stitching density on the infiltration pattern and filling time. The stitch row 
spacings were 6.35 mm, 12.7 mm, and 25.4 mm, respectively. 

The predicted flow front locations for the three preforms are shown in Figures 6.13 - 
6.15. In the figures, the color bands represent the flow front location at different times, and 
the units are in seconds. During the infusion process, the resin first moves across the high- 
permeable distribution medium and infiltrates the top carbon face sheets, and then ‘leaks’ 
through the stitch paths to wet out the bottom face sheets. 



Min 0 @Nd 608 


Figure 6.13: Predicted flow front progression in the stitched foam core preform with 6.35 
mm stitch row spacing 
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Max 1360 @Nd 5555 
Min 0 @Nd 637 

Figure 6.14: Predicted flow front progression in the stitched foam core preform with 12.7 
mm stitch row spacing 
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Max 1351 @Nd 3915 
Min 0 @Nd 523 

Figure 6.15: Predicted flow front progression in the stitched foam core preform with 25.4 
mm stitch row spacing 
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Both the digital images captured in the experiment and the flow front fringe plot predicted 
by the simulation show that the infiltration pattern on the top surfaces of the preform are 
quite uniform for all the three panels with different stitch row spacings. The predicted 
infiltration pattern in the top surfaces of the three preforms are shown in Figures 6.16 - 
6.18. Figures 6.19 through 6.26 show the experimentally observed progression of the flow 
front at the surface of the top face sheet of the preform with stitch row spacing of 12.7 
mm. Note that the flow is much slower along the edges of the panels because there is no 
distribution medium. Shown in Figures 6.27 - 6.29 are comparisons between the measured 
and predicted flow front locations as a function of time for the top surfaces of the preforms. 
The model prediction agrees with the experimental results fairly well. 
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Figure 6.16: Predicted flow front progression in the top surface of the stitched foam core 
preform with 6.35 mm stitch row spacing 
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Figure 6.17: Predicted flow front progression in the top surface of the stitched foam core 
preform with 12.7 mm stitch row spacing 
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Figure 6.18: Predicted flow front progression in the top surface of the stitched foam core 
preform with 25.4 mm stitch row spacing 



110 


130 



Figure 6.19: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 25.4 cm from injection edge 



Figure 6.20: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 30.5 cm from injection edge 
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Figure 6.21: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 35.6 cm from injection edge 







Figure 6.22: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 40.6 cm from injection edge 
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Figure 6.23: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 45.7 cm from injection edge 



Figure 6.24: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 55.9 cm from injection edge 
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Figure 6.25: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 61.0 cm from injection edge 



Figure 6.26: Top surface of stitched foam core preform with 12.7 mm stitch row spacing. 
The preform has been completely infiltrated 
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Figure 6.27: Measured and calculated flow along the top surface of the sandwich preform 
with 6.35 mm stitch row spacing 
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Figure 6.29: Measured and calculated flow along the top surface of the sandwich preform 
with 25.4 mm stitch row spacing 




The lower face sheet resin infiltration patterns were observed quite different for the pre- 
forms with different stitch spacings. Figures 6.30 through 6.36 show the flow front progression 
along the bottom surface of the panel with the stitch row spacing of 6.35 mm, Figures 6.37 
through 6.46 show the flow front progression along the bottom surface of the panel with the 
stitch row spacing of 12.7 mm, and Figures 6.47 through 6.51 show the flow front progres- 
sion along the bottom surface of the panel with the stitch row spacing of 25.4 mm. Resin 
enters the lower face sheet through the vertical stitches and then permeates normal to the 
stitching direction to fill the areas between the stitches. This results in the resin flow along 
the stitching leading the flow in the regions between the stitches. For the preform of 6.35 
mm stitch spacing, this lead-lag phenomenon is not evident and the flow pattern appears 
to be uniform due to the fact that the rows of stitches are close together and there is little 
lag in the area between the rows. When the row spacing increases, the lag length becomes 
larger and the infiltration pattern illustrates obvious nonuniformity. 



Figure 6.30: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 25.4 cm from injection edge 
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Figure 6.31: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 30.5 cm from injection edge 



Figure 6.32: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 35.6 cm from injection edge 
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Figure 6.33: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 40.6 cm from injection edge 



Figure 6.34: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 48.3 cm from injection edge 
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Figure 6.35: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
Flow front is approximately 61.0 cm from injection edge 



Figure 6.36: Bottom surface of stitched foam core preform with 6.35 mm stitch row spacing. 
The infiltration of the preform is nearly complete 
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Figure 6.37: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 25.4 cm from injection edge 



Figure 6.38: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 30.5 cm from injection edge 
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Figure 6.39: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 35.6 cm from injection edge 
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Figure 6.40: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 40.6 cm from injection edge 
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Figure 6.41: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 45.7 cm from injection edge 



Figure 6.42: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 50.8 cm from injection edge 


124 


144 



Figure 6.43: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 55.9 cm from injection edge 



Figure 6.44: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 61.0 cm from injection edge 
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Figure 6.45: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
Flow front is approximately 63.5 cm from injection edge 



Figure 6.46: Bottom surface of stitched foam core preform with 12.7 mm stitch row spacing. 
The preform has been completely infiltrated 


126 


146 





Figure 6.47: Bottom surface of stitched foam core preform with 25.4 mm stitch row spacing. 
Flow front is approximately 20.3 cm from injection edge 



Figure 6.48: Bottom surface of stitched foam core preform with 25.4 mm stitch row spacing. 
Flow front is approximately 28.0 cm from injection edge 
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Figure 6.49: Bottom surface of stitched foam core preform with 25.4 mm stitch row spacing. 
Flow front is approximately 35.6 cm from injection edge 



Figure 6.50: Bottom surface of stitched foam core preform with 25.4 mm stitch row spacing. 
Flow front is approximately 40.6 cm from injection edge 
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Figure 6.51: Bottom surface of stitched foam core preform with 25.4 mm stitch row spacing. 
Flow front is approximately 58.4 cm from injection edge 

Figures 6.52 - 6.54 show the model predicted infiltration patterns in the lower face sheets 
of the three preforms. The model accurately predicted the resin infiltration patterns observed 
in the experiments. When the stitch row spacing increases, the flow front shape in the bottom 
face sheet becomes more nonuniform. Extreme nonuniformity of the flow front shape could 
result in the formation of voids (dry spots) during infiltration process. This void entrapment 
occurs because there is more than one flow front, and when these fronts meet up, they will 
coalesce forming voids via entrapment. The possible void locations on the bottom face sheet 
of the sandwich preform with the stitch row spacing of 25.4 mm are predicted as shown in 
Fig 6.55. Fig. 6.55 illustrates the flow front shapes on the bottom surface of the preform 
at different time. Note that there are some spots near the ends of the preform, which have 
different color from the area around. This indicates that the resin fills these locations from 
around. Therefore, the resin fronts meet at these locations and form voids via entrapment. 

The measured and predicted total infiltration time for all the three preforms were com- 
pared in Table 6.3. Both the experiment and simulation show that the filling time increases 
with increasing stitch density. Initially one might assume that the preform with the highest 
stitch density will have the lowest infiltration time as a result of the large number of pen- 
etrations feeding resin to the lower face sheet. However, the panel with the highest stitch 
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Figure 6.52: Predicted flow front progression in the bottom surface of the stitched foam core 
preform with 6.35 mm stitch row spacing 
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Figure 6.53: Predicted flow front progression in the bottom surface of the stitched foam core 
preform with 12.7 mm stitch row spacing 
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Figure 6.54: Predicted flow front progression in the bottom surface of the stitched foam core 
preform with 25.4 mm stitch row spacing 
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Figure 6.55: Predicted flow front progression in the bottom surface of the stitched foam core 
preform with 25.4 mm stitch row spacing. Note the voids near the ends of the preform 
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Table 6.3: Measured and predicted infiltration time 


Stitch row 
spacing (mm) 

Infiltration time (min) 

Experiment 

measurements 

Model 

prediction 

% error 

6.35 

29.77 

25.59 

14.04 

12.7 

26.40 

22.67 

14.13 

25.4 

23.47 

22.52 

4.05 


density also requires a greater volume of resin to fill the face sheets and the needle penetra- 
tions in the foam core, which increases the overall infiltration time. The prediction for the 
infiltration times underpredicted the measured filling time by 4 to 14%. Differences might 
be attributed to the inaccurate permeability values used for the carbon face sheets. The 
permeabilities of the Tenax multiaxial warp knit preform were taken from data obtained in 
a previous investigation for different stitch densities. On the other hand, the presence of the 
foam core might change the permeabilities of the carbon fabric. An accurate permeability 
model will be the subject of the further investigation. 
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6.5 Conclusions 


Foam core sandwich preforms with through-the-thickness stitching were successfully resin 
infiltrated by the VARTM process. The through-the-thickness stitching permits resin infil- 
tration of the entire preform by placing a single injection tube along one edge of the upper 
face sheet. The resin first travels along the distribution medium and through the thickness 
of the upper face sheet. Then the resin travels through the foam core along the stitch pen- 
etrations and into the lower face sheet. The stitch density affects both the resin infiltration 
patterns and the total infiltration time. The infiltration front was uniform across the width 
of the upper face sheet. However, the flow patterns on the lower face sheet are nonuniform 
and clearly depend on the stitching spacing. The flow front along the stitching leads the 
flow front between the stitches. The lag length increases with increasing stitch row spacing 
and results in more nonuniform infiltration pattern in the bottom face sheet of the sandwich 
preform. Extreme nonuniformity of the flow front shape could result in the formation of the 
voids (dry spots). The total infiltration time decreases when the stitch spacing increases. 

The 3DINFIL computer code was used to simulate resin infiltration of the stitched foam 
core sandwich preforms. In the model, each row of stitches was described as a porous strip, in 
which the resin is only allowed to flow in through-the-thickness direction. The permeability 
of the strip was calculated analytically using a modified straight capillary model. Simulation 
results of the computer code were verified using flow visualization experiments. The model 
accurately predicted the resin infiltration pattern during VARTM of the stitched foam core 
preforms. The calculated total infiltration time underpedicted the measured filling time by 
4 to 14%. The differences may be attributed to the permeability values that are input into 
the model for the carbon face sheets. 
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Chapter 7 


Conclusions and Future Work 


7.1 Conclusions 

This investigation developed and verified a comprehensive Vacuum Assisted Resin Transfer 
Molding (VARTM) process simulation model. The model incorporates resin flow through 
the preform, compaction and relaxation of the preform, and viscosity and cure kinetics of the 
resin. The computer model can analyze the resin flow details, track the thickness change of 
the preform, predict the total infiltration time and final fiber volume fraction of the parts, and 
determine whether the resin could completely infiltrate and uniformly wet out the preform. 

Flow of resin through the preform is modeled as flow through porous media, Darcy’s 
law combined with the continuity equation for an incompressible Newtonian fluid forms the 
basis of the flow model. During the infiltration process, it is well accepted that the total 
pressure is shared by the resin pressure and the pressure supported by the fiber network. 
With the progression of the resin, the net pressure applied to the preform decreases as a 
result of increasing local resin pressure. This leads to the springback of the preform, and is 
called the springback mechanism. On the other side, the lubrication effect of the resin causes 
the rearrangement of the fiber network and an increase in the preform compaction. This 
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is called the wetting compaction mechanism. The thickness change of the preform is deter- 
mined by the relative magnitude of the springback and wetting deformation mechanisms. 
In the compaction model, the transverse equilibrium equation is used to calculate the net 
compaction pressure applied to the preform, and the compaction test results are fitted to 
give the compressive constitutive law of the preform. The Finite Element/Control Volume 
(FE/CV) method is adopted to find the flow front location and the fluid pressure. The code 
features the ability of simultaneous integration of 1-D, 2-D and 3-D element types in a single 
simulation, and thus enables efficient modeling of the flow in complex mold geometries. 

VARTM of two flat composite panels was conducted to verify the simulation model. The 
composite panels were fabricated using the SAERTEX multi-axial warp knit carbon fiber 
fabric and SI-ZG-5A epoxy resin. Panel 1 contained one stack of the carbon fabric, and 
Panel 2 contained four stacks of the fabric. The parameters verified included the flow front 
location and preform thickness change. For Panel 1, the flow front locations were accurately 
predicted while the predicted resin infiltration was much slower than measured for Panel 
2. The disagreement is attributed to the permeability model used in the simulation, which 
failed to consider the interface flow in the unstitched preform containing more than one stack 
of the fabric under very low compaction force. The predicted transverse displacements agree 
well with the experiment measurement qualitatively, but not quantitatively. The reasons 
for the differences were discussed, and further investigations are recommended to develop a 
more accurate compaction model. 

The simulation code was also used to investigate the VARTM of a new form of sandwich 
structure with through-the-thickness reinforcements, which is being considered for use in 
primary aircraft structure. The infiltration of three foam core sandwich preforms with dif- 
ferent stitch densities was studied. The objective of the study was to determine whether the 
preforms could be completely infiltrated and how the stitch density affects the infiltration 
process. The visualization experiments were conducted to verify the simulation. The model 
accurately predicted the resin infiltration patterns. The calculated filling times underpre- 
dicted experimental times by 4 to 14%. The model revealed the resin flow details and found 
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that increasing the stitch spacing shortens the total filling time, but increases the nonuni- 
formity of the flow front shape. Extreme nonuniformity of the flow front shape could result 
in the formation of the voids. 


7.2 Future Work 

Compaction and relaxation behavior of the preform was found to be very complex, nonlinear 
and loading history dependent. The compressive constitutive relation obtained by the com- 
paction test conducted in the quasi-static state was found not accurate enough to model the 
compaction behavior of the preform during the VARTM process. Further investigations are 
recommended to find a more accurate constitutive law representing the compaction response 
of the preform. Based on this investigation, a more accurate compaction model should be 
developed to predict the thickness change of the preform during infiltration process quanti- 
tatively. 
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1 Introduction 

This report presents the summary of the work perfonned at McGill University on the 
VARTM modeling and verification using 3DINFIL . The report is organized as follows: 

Section 2 presents the cure kinetic and viscosity models for the new VARTM resin system 
VR56-19. Section 3 compares the simulation of the infiltration of a unidirectional preform 
and the results measured experimentally. Section 4 discusses the modeling of flat and 
complex shapes with 3DINFIL. A sensitivity analysis is also presented where critical 
parameters that affect the infiltration are discussed. Finally, Section 5 compares the 
infiltration of a hat stiffened preform with the simulation results. 
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2 Matrix Characterization 

In order to meet variety of structural and manufacturing requirement, many different resin 
systems exists. The system of interest here is a VARTM two-part epoxy resin system that is 
suitable for structural aircraft composite applications. The resin to be analysed is the VR56- 
19, a derivative of the VR56-4 resin manufactured by Applied Poleramic Inc. The VR56-4 
resin has already been fully characterized for application in the VARTM process and 
simulation were carried out using the process model COMPRO [1]. The same 
characterization of cure kinetic and viscosity is to be done with the VR56-19. 

2.1 Cure Kinetics Model 

The objective here is to find the cure kinetics using a combination of isothermal and dynamic 
differential scanning calorimeter (DSC). All test were performed on a TA Instruments DSC- 
Q100 differential scanning calorimeter. The total heat of reaction was measured from the 
dynamic scans at l°C/minute from room temperature up to 250°C. The isothermal tests were 
performed at temperatures ranging from 60°C to 160°C. The isothermal test procedure was to 
rapidly heat the cell to the desired temperature and maintain it for up to 8 hours (in most 
cases). The cell was then rapidly cooled down to room temperature so that a dynamic scan at 
l°C/minute can follow to measure the residual heat of reaction. The following tables show 
the dynamic & isothennal DSC tests executed as well as the results. 


Table 1 VR56-19 dynamic DSC test matrix. 


Test 

# 

Test File Name 

Heating 

Rate 

Total 

HR 

Peek 

Temp. 

Freezer 

Time 

Sample 

Weight 

°C / min. 

J / g 

°C 

hours 

mg 

1 

DynamicScan_29-04-2004 

1 

650.2 

126.14 

12 

16.5 

2 

DynamicScan_01 -05-2004 

1 

585.9 

130.26 

0 

6.55 

3 

DynamicScan-Frozen_01 -05-2004 

1 

686.5 

126.78 

72 

10.69 

4 

DynamicScan_02-05-2004 

1 

613.3 

131.54 

0 

6.63 

5 

MDSC_DynamicScan_25-05-2004 

1 

530.5 

128.23 

0 

7.38 

6 

KIM_08- 13-2003 

1 

618.2 

170.35 

48 

2.28 
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Table 2 VR56-19 isothermal DSC test matrix. 


Test 

# 

Test File Name 

Iso. 

Temp. 

Iso. 

Time 

Residual 

Heat 

Freezer 

Time 

Sample 

Weight 

°C 

hours 

J / g 

hours 

Mg 

1 

Isothermal60_25-05-2004 

60 

8 

352.7 

24 

7.67 

2 

Isothermal60-2_31 -05-2004 

60 

8 

394.4 

72 

11.34 

3 

Isothermal80_27-05-2004 

80 

8 

133.4 

0 

5.78 

4 

Isothermal80-2_04-06-2004 

80 

8 

129.0 

2 

7.50 

5 

Isothermal 100_28-05-2004 

100 

8 

72.0 

24 

7.08 

6 

Isothermal 1 00-2_06-06-2004 

100 

8 

74.1 

72 

7.56 

7 

Isothermal 1 20_29-05-2004 

120 

8 

37.6 

48 

6.74 

8 

Isothermal 120-2_07-06-2004 

120 

6 

35.85 

0 

8.30 

9 

Isothermal 1 40_30-05-2004 

140 

8 

6.67 

0 

9.19 

10 

Isothermal 1 40-2_08-06-2004 

140 

5 

6.07 

24 

8.78 

11 

Isothermal 160_09-06-2004 

160 

3 

0 

36 

7.82 


Raw data from the DSC experiments consists of measurements of heat flow with respect to 
either time or temperature. The apparatus software is used to compute via integrals the total 
heat of reaction, or total residual heat. Those values are also in Tables 1 & 2. In order to 
better understand the behaviour of the epoxy resin and to characterize it properly, the DSC 
data are imported in excel for further analysis. Using pre-established formulas and raw DSC 
data, the degree of cure and cure rate can be computed. Two runs of all experiments are 
performed to validate the results. 

To visualize the curing dynamic, two sets of graphs are draw: -1- degree of cure with respect 
to time, -2- cure rate with respect to degree of cure. As a very similar resin (VR56-4) is 
already characterized, it is a great interest to compare the graphical behaviour of this resin 
with the VR56-19. The following 12 figures show on the comparison graphs -1- and -2- for 
each tests perfonned. 
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0 50 100 150 200 250 

Temperature (°C) 

Figure 1 Degree of cure variation for dynamic DSC tests. 
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Figure 2 Cure rate for dynamic DSC tests. 
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Figure 3 Degree of cure variation for 60°C isothermal cure. 



Degree of Cure 


Figure 4 Rate of cure variation for 60°C isothermal cure. 
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Figure 5 Degree of cure variation for 80°C isothermal cure. 



Figure 6 Rate of cure variation for 80°C isothermal cure. 


- 6 - 


176 


Rate of Cure (/s) Degree of Cure 


Final Report - VARTM Process Modeling and Verification 



Figure 7 Degree of cure variation for 100°C isothermal cure. 
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Figure 8 Rate of cure variation for 100°C isothermal cure. 
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Figure 9 Degree of cure variation for 120°C isothermal cure. 



Figure 10 Rate of cure variation for 120°C isothermal cure. 
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Figure 11 Degree of cure variation for 140°C isothermal cure. 



Figure 12 Rate of cure variation for 140°C isothermal cure. 
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Based on a visual analysis of the above figures, it can easily be seen that the VR56-19 
behaves quite differently than the VR56-4. This difference is mainly due to the fact that 
VR56-19 does not seem to have the two peaks in the isothermal graphs -2-. This is good in 
the way that the mathematical model of the resin should be simpler to find. Following is the 
procedure that is followed to find it. 

In order to find as quickly as possible the matching model and to optimize the parameters of 
that model for the DCS data, the difference between experimental and calculated rate of cure 
are minimized using Excel equation solver. In order to get results as accurate as possible, the 
difference between experimental and calculated degree of cure will also be considered, han 
the Excel solver is asked to minimize the sum of the two differences above mentioned. 
Curves from model 6 and 7, which are the ones that fit the data the most accurately, are 
compared with the experimental curves. It turns out that the model 6 yields better results, 
being based on a modified auto-catalytic equation accounting for a shift form kinetics to 
diffusion control. 


da _ Ka m (l - a)" 

dt l + e C {“"(«co+«crr)} 

K = Ae ( - AE/RT) 


( 1 ) 


The significance of the various terms in the equation used is presented in Table 3. The 
activation energy, AE, is calculated from the slope of the natural logarithm of the isothennal 
cure rate, In(daJdi), vs \/T at a number of different resin degrees of cure. All other model 
constants were detennined using Excel equation solver, using data from both isothermal and 
dynamic DSC measurements. To get the model to fit even better the results, the activation 
energy was also allowed to vary slightly when solving with Excel equation solver. 

Table 3 - Parameters for VR56-19 cure kinetics model. 


Parameter 

Value 

Activation energy 

AE = 54.95 kJ/gmole 

Pre-exponential cure rate coefficient 

A = 1.601x10 4 Is 

First exponential constant 

m = 0.4055 

Second exponential constant 

« = 1.362 

Diffusion constant 

C = 26.93 

Critical degree of cure at T=0°K 

&co = "0.2258 

Constant accounting for increase in critical resin 
degree of cure with temperature 

a CT = 2.762x1 0 -3 /K 


The following figures show the best model found up to now to fit the experimental data. It 
can be seen that at the model fits almost perfectly for temperatures of 80 °C to 120°C and 
differs slightly for temperatures of 60 °C or 140 °C. 
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Figure 13 Comparison between predicted and measured cure rate for a l°C/min ramp. 



Figure 14 Comparison between predicted and measured degree of cure for a l°C/min 

ramp. 
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Figure 15 Comparison between predicted and measured cure rate for isothermal cure. 



Figure 16 Comparison between predicted and measured degree of cure for isothermal 

cure. 
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2.2 Viscosity Model 

The VR56-19 viscosity rheometer results were already kn own at the beginning of this 
characterization research. They come from dynamic tests and a series of isothermal cure tests 
ranging from 60°C to 140°C. In those tests, the samples were sheared between two parallel 
discs until the resin reached its gel point. Again, in order to compare VR56-19 viscosity 
results with the VR56-4 ones, graphs of viscosity with respect to time are shown in figures 
17 and 18. It can be seen from those graphs that the viscosity of the VR56-19 is always a 
higher than the one of the VR56-4, but the overall behaviour if very similar. 
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Figure 18 Comparison between VR56-4 and VR56-19 isothermal viscosity tests. 

The viscosity model used in this characterization study is the same as the one of the VR56-4 
and is as follow: 


H = A h exp(£ /; / RT \a g /(a g - a^ A+Ba) (2) 

where A M , E M , A and B are experimentally determined parameters, R is the universal gas 
constant and a g is the degree of cure at gelation. The slope of a linear regression through the 
data of In /u versus 1 IT at low resin degree of cure ( a ~ 0) yields the value of E u . The data 
from the dynamic runs was used in this calculation. The gel point degree of cure ( a g ) was 
evaluated from the crossover point between the storage and the loss modulus (G’ and G”). 
Then, to evaluate the other constants (Table 4), a best fit was done using Excel equation 
solver with the constants A M , A and B to fit the experimental data from the isothermal and 
dynamic tests. 


Table 4 Parameters for VR56-19 viscosity model. 


Parameter 

Value 

Activation energy 

E m = 70.4 kJ/gmole 

Pre-exponential cure rate coefficient 

A m = 1.6x10" 12 Pa«s 

First exponential constant 

4 = 2.75 

Second exponential constant 

5=1.75 

Critical degree of cure at T=0°K 

a g = 0.54 
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This model capture the ambient viscosity, the reduction in viscosity due to temperature 
increase, and the rapid increase in viscosity as the resin reach the gel point. As it can be seen 
on the Figures 19 & 20, the model captures very well the curing viscosity in curing condition. 
The model does not capture very well the viscosity for the isothermal runs since there are no 
parameters in the model to set a minimum viscosity. The viscosity is then directly dependant 
on the temperature, as it can be seen in Figure 21. Even with those discrepancies on the 
isothermal runs, the model well captures the large increase in viscosity when getting close to 
the gel point. 
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Figure 19 Measured and predicted viscosity for a typical cure cycle. 
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Figure 20 Comparison between predicted and measured viscosity for dynamic cure 

condition. 
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Figure 21 Comparison between predicted and measured viscosity for isothermal cure 

condition. 
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3 Uniaxial Preform VARTM Modeling 

This section presents the results from simulations conducted using 3DINFIL for prediction of 
the infiltration of a preform made of the IM7 tackified carbon fiber uniaxial plies. The 
preform properties used for the simulations are derived from the permeability and 
compaction tests performed at NASA Langley. A series of models were generated following 
a run matrix and the results were compared to the results from a VARTM trial performed at 
NASA Langley (PANEL 012904). 

3.1 Preform characterization 

The perform characterization tests were conducted at NASA Langley using a special jig to 
measure in-plane and transverse permeability. An instrumented VARTM jig was used to 
measure the perform compaction behavior. 

3.1.1 Permeability 

Figure 22 shows the permeability data from the experiments conducted at NASA Langley. 
An eight plies layup ([-45/45/0/90] 2s ) tackified uniaxial carbon fibers was used for the tests. 
The permeability model in Equation 3 was fitted in the data and used in 3DINFIL 
simulations. 

IE-10 
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Figure 22 Permeability of [-45/45/0/90] 2s ) tackified uniaxial carbon preform. 
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K,, = 1.70xl0~ 13 (F f V 6 ' 3643 

/ \ (3) 

=1.80x10 15 (^) -3 - 5548 

Where K\\, K , and Vf are the in-plane permeability, transverse permeability and fiber 
volume fraction respectively. 

3.1.2 Compaction 

Figure 23 shows the compaction data from the experiments conducted at NASA Langley. 
An eight plies layup ([-45/45/0/90] 2 s) tackified uniaxial carbon fibers was used for the tests. 
The compaction model in Equation 4 was fitted in the data and used in 3DINFIL simulations. 
The initial porosity for the prefonn is 0.6829. 



♦ Dry - Experiment — 

— Dry - Model 

▲ Wet - Experiment ■ ■ 

■ ■ ■ Wet - Model 


Figure 23 Compaction of [-45/45/0/90] 2s ) tackified uniaxial carbon preform. 


£ dry 
wet 


0.4(l - e - 0M5p ) 


0.37 + 


0.09P 
16 + P 


(4) 


where Sdry, s wet and P are the dry preform strain, wet preform strain and the fiber pressure. 
The preform fiber volume fraction is computed internally by 3DINFIL using the following 
relation: 
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where (f) and s are the initial preform porosity and the preform strain {sd, y or £ wet ). 

Both penneability and compaction models were implemented in 3DINFIL as material 
IM 76KUNIT ACK_ 1 (ID = 16) in the database. 

3.2 VARTM Trial 

A VARTM infiltration of a 30 x 60 cm [-45/45/0/90] 2s tackified uniaxial carbon preform was 
performed at NASA Langley using a synthetic oil with a viscosity of 0.376 Pa.s (PANEL 
012904). The flow front position with time was measured on top and tool side and the results 
are presented in Figure 24. Due to the very small transverse permeability of the preform, an 
important time lag (200-300 seconds) between the top and tool side was observed. 
Furthermore, it was very difficult to measure the tool side flow front evolution due to the 
uneven flow front observed on the tool side. 



Figure 24 Flow front position during VARTM infiltration of [-45/45/0/90] 2s ) tackified 
uniaxial carbon preform with synthetic oil (viscosity of 0.376 Pa.s). 

3.3 Simulations 

Simulations of the infiltration of PANEL 012904 were perfonned to compare the flow front 
predicted by 3DINFIL and the experimental results shown in Figure 24. The finite element 
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mesh was generated using MSC PATRAN and the simulations were performed using 
3DINFIL Version 5.0 on a PC with the Linux operating system. 

3.3.1 Run definitions 


A 2-D mesh was chosen as the flow front was reasonably unifonn across the width of the 
panel. Figure 25 shows the geometric definition of the mesh with the different material 
regions and the location of the injection point. The distribution medium has a constant 
dimension with a length of 575 mm and a thickness of 1.5 mm. The preform length is 
constant (600 mm), but the thickness is varied in order to investigate the thickness effect on 
infiltration behavior. A constant number of 3500 quadrilateral four noded elements were 
used for all simulations. Typical simulation time was less than five minutes. 


Resin 



Figure 25 PANEL 012904 mesh definition. 

A constant injection pressure of 101000 kPa was applied at the top left node as shown in 
Figure 25 while all other nodes add an initial pressure set to zero. The material properties are 
summarized in Table 5. The medium and prefonn initial porosity was set to 75% and 
68.29% respectively. The resin viscosity was set to 0.376 Pa.s and was constant during the 
infiltration. The simulations used the flow and compaction solution algorithm. 

Table 5 Material properties. 



Medium 

Preform 

Initial porosity 

75% 

68.29% 

K xx 

8.3e-9 


Kyy 

3.0e-l 1 

Equation 3 

K zz 

8.3e-9 


Compaction 

/ 

Equation 4 


A total of seven runs were defined as shown in Table 6. The prefonn thickness was varied 
from 4.58 to 1.8 mm. The range was chosen from the thickness measured during dry and wet 
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compaction tests. A thickness of 4.58 mm corresponds to the initial preform thickness 
measured for the dry compaction tests. The minimum dry preform thickness measured under 
full vacuum was 2.75 mm, while the minimum thickness for a wet preform under full 
vacuum was measured at 2.53 mm. We also decided to investigate preform thickness below 
the minimum measured thickness to account for is possible errors in the prefonn thickness 
measurements 

A typical mesh is shown in Figure 26 It is clear that the preform is very thin compared to the 
length; therefore the mesh density was chosen to maintain an adequate aspect ration for the 
elements. Two elements were used across the medium thickness, while six elements were 
used across the prefonn thickness. 


Table 6 Runs matrix. 


Run Name 

Preform Thickness, t (mm) 

Run 1 

4.58 

Run 2 

4.0 

Run 3 

3.5 

Run 4 

3.0 

Run 5 

2.5 

Run 6 

2.0 

Run 7 

1.8 
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Figure 26 Typical mesh showing element aspect ratio used. 

3.3.2 Results 

Figure 27 shows the results for the predicted flow front evolution when the initial preform 
thickness is set to 4.58 mm (Run 1). Good agreement with the experiment was observed for 
the top side. However, a huge difference (1000 seconds) was observed for the tool side. 
Figure 28 and Figure 29 show the flow front predicted for all runs for the top side and tool 
side respectively. From these figures, it appears that the prefonn with a thickness of 2 mm 
(Run 6) matches the experiment. In this case, difference between the experiments was less 
than 10% on average. Figure 30 clearly illustrates the thickness effect on the infiltration 
time. Thickness has a string effect on the flow at the tool side, but surprisingly has also a 
slight effect on the flow on the top side. We believe that the interaction between the in-plane 
and through thickness flow is responsible for that behavior. For a thin panel, the flow gets to 
the bottom faster. Therefore, the advance of the flow in the preform, particularly at the top, is 
faster than a thicker preform. 
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Figure 27 Comparison between predicted and measured flow front position with a 
preform initial thickness of 4.58mm (Run 1). 
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Figure 28 Comparison between predicted and measured flow front position on the top 

side. 
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Figure 29 Comparison between predicted and measured flow front position on the tool 

side. 
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Figure 30 Effect of preform thickness on flow front time at the end of the medium 

position = 575mm). 
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3.4 Discussion 

The panel investigated in this work is very challenging to simulate because of the large 
difference between the in-plane and through the thickness permeability (three orders of 
magnitude). 3DINFIL simulations clarify the different phases during a VARTM infiltration: 
the filling of the medium (Phase I), the end of the medium filling (Phase II) and the end of 
the preform filling (Phase III). In Phase I, we observe that the medium is quickly infiltrated 
with some resin flowing in the preform. But, as soon as the medium is full (Phase II), the 
resin pressure increases very rapidly and the flow rate increases in the preform. In Phase III, 
the flow front progresses slowly as the end of the preform is infiltrated in-plane. 

Two phenomena are in competition: decreasing the medium length to accelerate the flow in 
the preform (pressure effect) also increases the length of the medium free preform where 
flow rates are low (Phase III). VARTM simulations could help finding the optimum medium 
length to ensure complete prefonn infiltration in the processing window defines by the resin 
curing behavior. 

3.5 Conclusions 

In this study, we used 3DINFIL to simulate the resin flow through a uniaxial tackified carbon 
fiber preform. The compaction and penneability models derived from experiments on the 
same preform were used in the simulations. We mainly found that: 

• Infiltration time is very sensitive to prefonn thickness 

• Current implementation of compaction-flow does not account for thickness variations 
during the infiltration process 
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4 3DINFIL Modeling 

4. 1 Preliminary sensitivity analysis 

In this first study, a 2D simple geometry was considered. The objective of this study was to 
stress the importance of accurate physical properties measurements. The panel modeled is 
600mm long and 4.58mm thick (Figure 31). However the thickness is changing during the 
infusion. Above the fiber prefonn, a medium layer of 575mm by 1.5mm is set. Finally, all 
the perform-medium assembly is embedded with a plastic bag not considered here. 


Resin 



Figure 31 Model geometric definitions. 

Quadrilateral finite elements (Figure 32) were used to mesh the geometry. A total number of 
3500 elements were used in the simulations. This number enables to get accurate results 
while keeping a calculation time below 5 minutes. 



Figure 32 Typical mesh. 
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4.1.1 Simulations definition 

As for all the other parameters, the material properties were obtained from the measurements 
on the material used. The nominal material properties are given in Table 7. 


Table 7 Nominal material properties. 



Medium 

Preform 

Initial porosity 

75% 

68.29% 

K xx (m 2 ) 

8.3e-9 

2.543e-10 

Initial Permeability K yy (m 2 ) 

3.0e-ll 

1.068e-13 

K zz (m 2 ) 

8.3e-9 

2.543e-10 


The prefonn has been defined as a proper material with the permeability law defined by 
Equation 3. 

A constant resin inlet pressure at the top left of the mesh was set and a vacuum at all the 
other nodes is applied. In 3DINFIL, this is accomplished by setting a heat source of -1 (the 
negative signs refers to an inlet resin group, and the number only refers to the number of the 
group) and a heat flux of 1 (initial fill factor) at the first node, and a heat source of 0 (initial 
fill factor) in all the other nodes. 

The geometrical model was set as given in Figure 31. The medium and prefonn porosity 
were set to 68.29% and 75% respectively. 

In the general 3-D case, the permeability tensor has 6 independent components (3 by 3 
symmetric matrix). In our simple example here, this matrix is reduced to 2 constants for each 
material. The penneability used is defined in Table 8. 
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Table 8 Runs definition and material permeabilities. 



Medium 

Preform 

Run number 

Kxx (m 2 ) 

Kyy (m 2 ) 

Kxx (m 2 ) 

Kyy (m 2 ) 

1 (Nominal) 

8.300E-09 

3.000E-11 

2.543E-10 

1.068E-13 

2 

5.000E-09 

3.000E-11 

2.543E-10 

1.068E-13 

3 

1.000E-09 

3.000E-11 

2.543E-10 

1.068E-13 

4 

1.600E-08 

3.000E-11 

2.543E-10 

1.068E-13 

5 

8.300E-09 

1.500E-1 1 

2.543E-10 

1.068E-13 

6 

8.300E-09 

6.000E-1 1 

2.543E-10 

1.068E-13 

7 

8.300E-09 

9.000E-1 1 

2.543E-10 

1.068E-13 

8 

8.300E-09 

8.000E-12 

2.543E-10 

1.068E-13 

9 

8.300E-09 

3.000E-11 

5.086E-10 

1.068E-13 

10 

8.300E-09 

3.000E-11 

2.543E-09 

1.068E-13 

11 

8.300E-09 

3.000E-11 

1.272E-09 

1.068E-13 

12 

8.300E-09 

3.000E-11 

2.543E-1 1 

1.068E-13 

13 

8.300E-09 

3.000E-11 

2.543E-10 

2.136E-13 

14 

8.300E-09 

3.000E-11 

2.543E-10 

1.068E-12 

15 

8.300E-09 

3.000E-11 

2.543E-10 

5.340E-14 

16 

8.300E-09 

3.000E-11 

2.543E-10 

1.068E-14 


The main idea of this study is to vary one parameter while keeping the others constant to 
their nominal value. The reference point is however different for runs 1 to 8 and 9 to 16. The 
preform characteristics are kept constant in runs 9 to 16. 

4.1.2 Results 

Figure 33-36 show the results of this penneability study. 
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Figure 33 Influence of the medium permeability for the top side of the panel. 
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Figure 34 Influence of the medium permeability for the tool side of the panel. 
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Figure 35 Influence of the preform permeability for the top side of the panel. 
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Figure 36 Influence of the preform permeability for the tool side of the panel. 

The shape of the curves obtained confirms the linearity of the equations and the convergence 
of 3DINFIL code. Nevertheless, 3DINFIL may diverge particularly for low penneabilities 
(Figure 36, Run 16). Actually, in Run 16, the convergence conditions were reached. The 
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results may radically change from one simulation to the other. For instance, from Figure 36, 
an order of magnitude change in penneability can double the flow time. 

In general, the measurement of permeability is more important for the preform than for the 
preform. Indeed, the curves are more dispersed when the preform permeability is varied. 
However it is important to keep in mind that the medium is much thinner than the preform 
(ratio length/thickness is 3 times higher for the medium). 

The longitudinal permeability affects the slope of the curves, while the transverse 
penneability changes the gap between the top and tool response. 

The medium response is very sensitive to variations in longitudinal permeability, while the 
preform response is more affected by the transverse penneability value. This underlines the 
physical principle of infusion: the low longitudinal prefonn penneability is compensed by the 
medium. 

4.1.3 Conclusion 

The modeling of resin flow requires several physical constants: 

- Porosity ( 1 constant) 

- Penneability tensor (up to 6 independent components) 

- Geometry (>2 components for the simplest 2D model) 

This gives a minimum number of 4 physical properties. And this number quickly gets 
increasingly high, as more complex or realistic models are considered. 

The porosity is simple to evaluate. In this study, the importance of only two components of 
the permeability tensor was investigated. 

The two points to remember form this study are that having a representative model requires 
accurate measurements of permeability. The second one is that the essence of the resin 
transfer molding is to compensate low longitudinal penneability by adding a highly 
permeable medium. 

As we have seen, there are still many variables that rule the physics of a resin flow. Other 
parameters, like geometry, are as important to get accurate simulation results. 

4.2 Complex geometries simulations 

Simple 2D geometries are useful to understand the general concepts of infusion. However, 
the goal of processing a simulation is to predict the infusion of real parts. In this section, two 
examples of complex composite structures are presented. As the equations governing resin 
flow through a solid preform are complex, solving the problem with finite elements can be 
easily accomplished. In this section, the model of an angle laminate is presented. 
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This study is related to the investigation of the dimensional stability of a C-shaped composite 
material. When an angle piece part is processed, angular directions arise when the part is 
removed from the mold. The goal is to investigate the relationship between the infusion 
scheme and the part distortion. 

4.2.1 Model description 

This study is based on a C-shaped part presented in [5]. The general geometrical data is 
represented in Figure 37. The part is composed of a 225mm long preform and a 210mm 
medium. It is designed on a 100mm large square mold. The angles of the mold have a radius 
of 6.25mm. Finally the preform thickness is 3.6mm and the medium is 1.5mm thick. The 
preform is centered on the tool. The medium may be set in a centered or asymmetric 
configuration as shown in Figure 38. 
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Figure 37 C-shape geometry. 
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Figure 38 Influence of the preform permeability for the top side of the panel. 

4.2.2 Material properties 

The definition of the material properties is the most difficult task for the curved geometry, 
since the change of orientation has to be considered. Different sets of properties for each 
material orientation along the curve have to be defined. The basic physical properties of the 
medium and the fiber are given in the Table 9: 

Table 9 Material properties of the C-shaped. 



Medium 

Fiber 

Designation 

User defined 

Multiaxial Wrap Knit, 
unstitched, Tenax 1 stack 

Porosity 

0.75 

0.5 

Kxx (m 2 ) 

8.300E-09 

3.934E-1 1 

Permeability Kyy (m 2 ) 

tensor 

3.000E-1 1 

1.646E-1 1 



Kzz (m 2 ) 

8.300E-09 

9.418E-13 


For each orientation, the rotation tensor was defined to get the right axial and transverse 
properties following the relation below: 


K 

= curved 


cos 6 sin 0 0 

-sin# cos 6 0 K 

= straight 

0 0 1 

6-rotation tensor 


( 6 ) 
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where 6 represents the orientation angle of the element in Figure 39. 

This limits the complexity of the model. Indeed, four linear segments with four orientations 
approximated the corner. For two angles and two materials, this means 16 different material 
properties (of 12 data number each) had to be defined only for the curves. This is a 
significant limitation of 3DINFIL for curved models. However, the present 4-orientations’ 
approximation enabled to obtain interesting results. 



Figure 39 Meshing properties of the angle part 

4.2.3 Mesh and boundary conditions 

Contrary to the material properties case, to mesh a curved geometry is not so harder than the 
straight case; Patran’s meshing algorithm does all the work. A quadrilateral type of elements 
was chosen. To keep short calculation time, and as we were already limited by the number of 
elements in the curved parts, the number of elements was limited to 1200 (Figure 40). This 
simplification does not affect the accuracy of the results. 
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In order to study two configurations (centered and asymmetric), two sets of boundary 
conditions were modeled. However the only thing changing is the position of the resin inlet 
(Figure 38). The principle of the boundary conditions remains the same. 

The filling factors are not really different from the previous simulations. Nodes at resin inlet 
are set to a negative temperature of -1 and a heat source of 1 (filling factor). In all the other 
nodes a null heat source (no initial filling) to model the vacuum was used. 

4.2.4 Asymmetric infusion results 

The first model was the asymmetric configuration. It enables the observation of the main 
phenomena as the flow process or the influence of the two corners compared to a straight 
configuration. The time flow along the part is shown in Figure 41. The form of the response 
is quite identical to the straight model: the resin flows quickly along the medium and slows 
down as it reaches the end. The same diffusion gap through the thickness, although in this 
case, the part is much thinner and permeable compared to the panel used in Section 3.1. 



Curvilinear coordinate (cm) 


Figure 41 Time flow for the external and tool side of the angle part. 

One of the other similarities with the flat panel case is the succession of flow phases (Figure 
42). The medium fills rapidly (Phase I), then the flow intensifies (Phase II) and finally the 
flow slows down in the final part of the preform (Phase III). This sequence is rather clear in 
the infusion animation in Patran. This phenomenon can be seen by plotting the pressure 
evolution at the center of the fibers (Figure 42). 
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Time (s) 


Tool side 

External side 


Figure 42 Pressure for the center top point of the fiber (external and tool side): Three 
different phases of filling. 

As the simulations gave interesting results with a short calculation time, a 3D model was 
created to study the edge effects. This model only lied in an extrusion in the z-axis of the 
2D's one. However, as can be seen in Figure 43, the edge effects are negligible The flow is 
constant in the z-coordinate. The 3D case is only an extrusion of the 2D flow. Therefore it 
was not worth to consider 3D complexity and the 2D representation was used from this point. 
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Figure 43 3D model's resin pressure results: The code does not show important 

transverse effect. 
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4.2.5 Effect of the infiltration location 

The centered configuration gives interesting results about the infiltration time. Indeed, the 
simulations show that the part is fully infiltrated in about 30 seconds while it requires 1 
minute and 45 seconds with the asymmetric process (same order of magnitude as the 
experiments). From a production point of view, it seems to be more interesting to use the 
centered model. However when the fiber volume fraction is important, the choice is not as 
trivial (Figure 44). 



Curvilinear coordinate 


Asymmetric 

(Tool side) 

Asymmetric 

(External side) 

Centered (Tool 

side) 

Centered 

(External side) 


Figure 44 Fiber fraction along the piece for the two resin inlet configurations. 

The fiber volume fraction variation through the thickness is negligible in this model as seen 
in Figure 44. The two inlet cases show the same distribution of fiber fraction between the 
resin inlet and the medium end. The quantity of resin is more important closer to the inlet 
point with a gradient greater for the centered configuration that could lead to heterogeneous 
shape near the end. Thus it is difficult to decide what inlet technique to use. Only a residual 
stress analysis could help in this comparison. 

4.2.6 Effect of the medium length 

Finally a series of simulations were conducted to study the influence of the medium length 
for the asymmetric case. The results for the time flow are shown in Figure 45. 
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18.4125 cm 

19.2750 cm 

20.1375 cm 
21.0000 cm 

21.8475 cm 

22.6950 cm 


Figure 45 Time flow for different medium lengths (measure on the tool side). 

Contrary to some other studies, in this case the time response is monotone towards the 
medium length: the longer the medium is, the quicker the piece is infused. Here the fiber 
permeability is too high to see a real optimization of the medium length (in some infusion 
cases, the longest medium does not give the lowest infusion time). The order of magnitude of 
the resin velocity is rather the same through the liber and through the medium. So the piece 
will always get infused to the tool side. 
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Figure 46 Fiber fraction repartition at the final time of infusion for different medium 

lengths. 

The results for the fiber fraction confirm the previous results. As the medium is longer, less 
resin is required to fill the whole part and the fiber volume fraction is more homogeneous. 
This explains the offset noticed in Figure 46. 

4.2.7 Conclusion 

This study shows many interesting results concerning the simulation of the infusion process 
for curved geometry. Although the results are for an angle part they could be generalized to 
any geometry. First of all, the complexity of the material properties definition for each 
orientation is an important limitation of the 3DINFIL code. 

Nevertheless, interesting observations were made and order of magnitude for he solution 
could be achieved. Additional experiment data would be necessary to be able to compare the 
simulation results. 

The next step now should be to use the output of this model and to run a residual stress 
analysis to explain the angle variation phenomenon. 
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4.3 3DINFIL Sensibility analysis 

Computer simulations of manufacturing processes are increasingly used since it can save 
time and money. For VARTM modeling, it is very important to input accurate properties and 
boundary conditions in order to obtain valid results form simulation software. The objective 
of this study is to evaluate the effect of variation in material and processing parameters on the 
total infusion time of a flat panel. The results of this research will highlight the most 
significant parameters involved in the VARTM manufacturing process. 

4.3.1 Model description 

The model used in 3DINFIL VARTM software is a simple flat panel. The general process to 
manufacture such a panel is shown in the following figure: 


Spiral Tubing Distribution Media 



Figure 47 VARTM manufacturing process. 

Assuming that the resin flow starts unifonnly from the inlet spiral tubing, proper 
approximation is to model only a longitudinal cross-section of the panel (Figure 48). 
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Figure 48 Longitudinal cross-section of a VARTM panel. 

The resin inlet is assumed as a simple node at the left hand side of the medium. The 
distribution medium does not extend completely to the end of the panel, in order to reduce 
the flow front lag between the top and bottom of the preform. Figure 49 clearly shows the 
mesh geometry generated in Patran. 



Figure 49 - Geometry mesh for the VARTM panel model. 
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4.3.2 Sensitivity analysis parameters 

In this study, nine different parameters are varied by ±10% of the nominal values obtained 
from material characterization or actual processing conditions. The parameters are classified 
in two categories: material properties and processing condition (see below). 

Material Properties: 

In-Plane Preform Penneability (A) 

Through Thickness Preform Permeability (B) 

In-Plane Medium Permeability (C) 

Through Thickness Medium Permeability (D) 

Dry & Wet Prefonn Compaction (E) 

Resin Viscosity (F) 

Processing Conditions: 

Pressure Gradient (G) 

Medium Gap Length (H) 

Preform Thickness (I) 


Table 10 lists the input of the various parameters for the nominal case, the lower bound (- 
10%) and the upper bound (+10%). The VARTM processing variable investigated is the 
total infusion time, which corresponds to the time taken by the resin to reach the lower (tool 
side) right corner of the preform. The nominal values for all the parameters resulted in a good 
agreement between experimental result and the 3DINFIL simulation (Figure 50). 



0.00 0.10 0.20 0.30 0.40 0.50 0.60 

Flow Front Position (m) 

Figure 50 Comparison between experiment and 3DINFIL results for nominal case. 
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Table 10 Material and process related parameters. 



Constant 

Name 

Lower Bound 
Level 1 -10% 

Upper Bound 
Level 2 +10% 

Nominal 

(Ref.) 

A 

In-Plane Preform 
Permeability 

a II 

1.52E-13 

1.86E-13 

1.69E-13 


b || 

-6.3590 

-6.3590 

-6.3590 

B 

Through Thickness 
Preform Permeability 

a ^ 

1.64E-15 

2.00E-15 

1.82E-15 


b i 

-3.5544 

-3.5544 

-3.5544 

C 

In-Plane Medium 
Permeability 

... 

7.47E-09 

9.13E-09 

8.30E-09 

D 

Through Thickness 
Medium Permeability 

... 

2.70E-1 1 

3.30E-1 1 

3.00E-1 1 



constl dn , 

0.3732 

0.4251 

0.39913 



const 2 dn , 

-0.0436 

-0.0497 

-0.04667 

E 

Dry & Wet Preform 
Compaction 

constl„ et 

0.3324 

0.3786 

0.35551 


const 2 „ et 

0.0906 

0.1032 

0.09687 



const3 we , 

5.9807 

6.8122 

6.39643 




0.66587 

0.70001 

0.68294 

F 

Resin Viscosity 
(Pa-s) 

... 

0.3384 

0.4136 

0.376 

G 

Pressure Gradient (Pa) 

— 

90900 

111100 

101000 

H 

Medium Gap Length (m) 

— 

0.0225 

0.0275 

0.025 

■ 

Preform Thickness (mm) 

— 

1.80 

2.20 

2.00 


4.3.3 Description of material related parameters 


Preform Permeability (A & B): The preform permeability is calculated in 3DINFIL with 

Equation 7: 
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*i=«n cf 

K ± =a ± (v r ) b 


(7) 


where K\\, K , , and F/ are the in-plane, through thickness permeability and fiber volume 
fraction respectively. Constants a», a L , and b ± are obtained from permeability tests. 

Figure 51 shows the preform permeability range studied (Table 10) where the upper and 
lower bound correspond to variability observed in the permeability measurements. 



Fiber Volume Fraction 
Figure 51 Preform permeability range. 


Medium Permeability (C & D): The medium permeability is described by a constant in 

3DINFIL . Through thickness permeability includes the effect of the peel ply layer. 


Preform Compaction (E): The preform compaction is calculated in 3DINFIL with 

Equation 8: 
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|| _ ^const2 dry P \ 


s d ,y =const\ dry \\-e 
s = cons , I + . const2 - P 


const! wet + P 


( 8 ) 


where £d, y , £ wet and P are the dry preform strain, wet preform strain and the fiber pressure. 
Constants constld, y , const2d, y , constl wet , const2 wet and const3 wet are obtained from preform 
compaction tests. The preform fiber volume fraction is computed internally by 3DINFIL 
using the following relation: 



1 zt 

1 - £ 


( 9 ) 


where (j) and £ are the initial preform porosity and the preform strain (£d, y or £ wet ). It is 
important to note that the initial prefonn porosity was also changed for the lower and upper 
bound (Table 10). Figure 52 shows the range in dry and wet compaction curves defined by 
the lower and upper bounds. 
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Figure 52 Preform compaction range. 
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Resin Viscosity (F): The resin viscosity corresponds to a constant in 3DINFIL and is 

assumed constant throughout the entire infusion time (which may not be the case in a real 
experiment if the infusion time is long compared to gel time). 

4.3.4 Description of processing condition parameters 


Pressure Gradient (G): This parameter corresponds to the pressure differential between the 
inlet (resin) and outlet (vacuum). 


Medium Gap Length (H): This parameter is shown by the “Gap” dimension in Figure 48. 


Preform Thickness (I): This parameter is shown by the “Thickness” dimension in Figure 

48. 


4.3.5 Sensitivity analysis results 

In order to optimize the sensitivity analysis for nine parameters, a Taguchi design of 
experiments approach was selected. Since this is a 9 two-level factors experiment study, a 
detailed analysis can be performed using only 12 experiments (L-12 array). The Taguchi L- 
12 orthogonal array is used to design experiments with 8, 9, 10 or 11 two-level factors. It is 
good to note that 512 (2 9 ) different experiments would have to be perfonned in the case no 
design of experiments techniques where used. Table 1 1 shows the experiments description as 
well as the total infusion time obtained in 3DINFIL. 

The results for each simulation are evaluated in the Qualitek-4 software. The software 
provides the results for the main effects, the optimum condition and performance and an 
analysis of variance. The main effects give the infusion time if only one parameter was 
changed at a time and all others were kept constant (Table 12). The optimum condition and 
performance, which is not of interest in this study, gives the level to choose for each 
parameter in order to get the lowest possible infusion time, as an example. Then, the analysis 
of variance (ANOVA) describes the steps to follow in order to get the percentage effects of 
each parameter on the total infusion time, as seen in Table 13. 
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Table 11 Taguchi L-12 orthogonal array with simulation results. 


Experiments 

A 

B 

c 

D 

E 

F 

G 

H 

1 

Infusion 
Time (sec.) 

Nominal Case 










485 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

540 

2 

1 

1 

1 

1 

1 

2 

2 

2 

2 

771 

3 

1 

1 

2 

2 

2 

1 

1 

1 

2 

432 

4 

1 

2 

1 

2 

2 

1 

2 

2 

1 

299 

5 

1 

2 

2 

1 

2 

2 

1 

2 

1 

412 

6 

1 

2 

2 

2 

1 

2 

2 

1 

2 

590 

7 

2 

1 

2 

2 

1 

1 

2 

2 

1 

455 

8 

2 

1 

2 

1 

2 

2 

2 

1 

1 

334 

9 

2 

1 

1 

2 

2 

2 

1 

2 

2 

583 

10 

2 

2 

2 

1 

1 

1 

1 

2 

2 

622 

11 

2 

2 

1 

2 

1 

2 

1 

1 

1 

568 

12 

2 

2 

1 

1 

2 

1 

2 

1 

2 

333 

Total 

18 

18 

18 

18 

18 

18 

18 

18 

18 



* 1 = Lower Bound (-10%) 
2 = Upper Bound (+10%) 
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Table 12 Main effects of each parameter on total infusion time from Taguchi analysis. 


Factors 

Infusion Time (s) 

Lower 

Bound 

Effect 

Higher 

Bound 

Effect 

In-Plane Preform Permeability 

507.4 

71 4.5% 

482.5 

^ 0.6% 

Through Thickness Preform 

Permeability 

519.2 

7! 7.0% 

470.6 

'J 3.0% 

In-Plane Medium Permeability 

515.6 

71 6.2% 

474.3 

5J 2.3% 

Through Thickness Medium 

Permeability 

502.2 

7! 3.5% 

487.7 

7! 0.4% 

Dry & Wet Preform Compaction 

591.1 

71 21.8% 

398.7 

il 17.9% 

Resin Viscosity (Pa-s) 

447.0 

'J 7.9% 

542.8 

7! 11.8% 

Pressure Gradiant (Pa) 

526.2 

71 8.4% 

463.6 

N 4.5% 

Medium Gap Length (m) 

466.2 

il 4.0% 

523.7 

7! 7.9% 

Preform Thickness (mm) 

434.6 

ii 10.5% 

555.2 

71 14.4% 
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Table 13 Analysis of variance of parameters on infusion time from Taguchi analysis. 


Factors 

DOF 

Sums of 
Squares 

Variance 

F-Ratio 

Pure Sum 

Percent 

In-Plane Preform 

Permeability (m 1 2 3 4 5 6 7 8 9 ) 

1 

1,860.0 

1,860.0 

3.5 

1,334.3 

0.6 

Through Thickness Preform 
Permeability (m 2 ) 

1 

7,076.2 

7,076.2 

13.5 

6,550.4 

3.0 

In-Plane Medium 

Permeability (m 2 ) 

1 

5,117.1 

5,117.1 

9.7 

4,591.3 

2.1 

Through Thickness Medium 
Permeability (m 2 ) 

1 

630.4 

630.4 

1.2 

104.7 

0.1 

Dry & Wet Preform 
Compaction 

1 

111,053.1 

111,053.1 

211.2 

110,527.3 

50.3 

Resin Viscosity (Pa-s) 

1 

27,552.1 

27,552.1 

52.4 

27,026.3 

12.3 

Pressure Gradiant (Pa) 

1 

11,768.6 

11,768.6 

22.4 

11,242.9 

5.1 

Medium Gap Length (m) 

1 

9,918.8 

9,918.8 

18.9 

9,393.0 

4.3 

Preform Thickness (mm) 

1 

43,609.0 

43,609.0 

82.9 

43,083.2 

19.6 


Error 2 1,051.4 525.7 2.6 


Total: 11 219,636.6 100.0 

From the analysis of variance, we can easily see that some parameters, such as perform 
compaction, are very important to understand before running computer simulation. In order 
of importance (Figure 53) are: 

1- Dry & Wet Preform Compaction (50.3 %) 

2- Preform Thickness (19.6 %) 

3- Resin Viscosity (12.3 %) 

4- Pressure Gradiant (5.1 %) 

5- Medium Gap Length (4.3 %) 

6- Through Thickness Prefonn Permeability (3.0 %) 

7- In-Plane Medium Permeability (2.1 %) 

8- In-Plane Preform Permeability (0.6 %) 

9- Through Thickness Medium Permeability (0.1 %) 
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Figure 53 Parameters overall effect on total infusion time. 

There is also an error of 2.64 % on the results, which is caused by the fact that the L-12 
orthogonal array used to solve the Taguchi system was not complete, as we used only 9 of 
the 1 1 possible columns. 

4.3.6 Conclusion: 

This case study first showed that a well-defined composite process model could provide 
accurate predictions of infusion time. This study also shows that it is very important to well 
understand the individual effect of each parameter used in 3DINFIL. It is clear that small 
variations in the prefonn compaction behaviour can have a significant effect on the preform 
infiltration time. Therefore, great care must take to measure the preform compaction 
behaviour and particularly the prefonn initial thickness. 
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5 Hat Stiffened Panel Simulations 

A simulation of the infiltration of a hat stiffener panel was performed to compare the flow 
front predicted by 3DINFIL and the experimental results. The finite element mesh was 
generated using MSC PATRAN and the simulations were performed using 3DINFIL Version 

5.0 on a PC with the Linux operating system. 

5.1 Run definitions 

A 2-D mesh was chosen to simplify the model. The flow front was reasonably uniform 
across the width of the panel but two sets of experimental results were considered. One in 
the middle of the panel and one in the last section to be infused (quicker infusion), as shown 
in Figure 54. 



Figure 54 Two sets of experimental results following black lines. 

Figure 55 shows the geometric definition of the hat stiffened panel with the different material 
regions, dimensions and the location of the injection point. The gap between the end of the 
MAWK perform and the end of the medium is initially set to 1 inch (25.4 mm). 


- 52 - 


222 



Final Report - VARTM Process Modeling and Verification 



_ 

/ 


Figure 55 - Hat stiffened panel geometry and materials. 

A constant injection pressure of 101,000 Pa was applied at the top left node as shown in 
Figure 55 while all other nodes add an initial pressure set to zero. The material thicknesses 
used in the model are summarized in Table 14 and the properties Table 15. Motor oil with a 
constant viscosity of 0.376 Pa.s was used for the experiment so the resin viscosity of the 
model was also kept constant at this same value during the infiltration. The simulations used 
the flow and compaction solution algorithm. 

Table 14 Material thicknesses used in the 3DINFIL model. 


Material 

Thickness (mm) 

Medium 

1.500 

1 Stack MAWK 

1.850 

1 Layer Biaxial 

0.318 

2 Layer Biaxial 

0.636 


Table 15 Material properties used the 3DINFIL model. 
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Medium 

MAWK 

Biaxial 

Initial Porosity 

0.752 

0.5714 

0.7507 

Permeability 

Model 

Constant 

None 

a 

b 

a 

b 

K xx 

8.30E-09 

2.00E-13 

-7.620 

6.21E-13 

-5.836 

Kyy 

5.00E-1 1 

6.00E-15 

-7.294 

6.21E-13 

-5.836 

K zz 

8.30E-09 

4.00E-13 

-5.363 

6.21E-13 

-5.836 

Compaction 

dry_const1 

None 

0.151 

0.475 

dry_const2 

-0.034 

0.065 

wet_const1 

None 

0.122 

0.306 

wet_const2 

0.035 

0.211 

wet_const3 

14.076 

3.168 


A typical mesh is shown in Figure 56. It is clear that the preform is very thin compared to 
the length; therefore the mesh density was chosen to maintain an adequate aspect ratio for the 
elements while keeping a low number of elements. 
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Figure 56 Typical mesh showing element aspect ratio used. 

5.2 Results 

Figure 57 shows the point used to compare the infusion time in 3DINFIL with the 
experimental results. Basically we consider the time the flow front takes to reach each of 
these points. Points 1-13 are associated with the medium side and points 14-22 are 
associated with the tool side. 



Figure 57 Points used for the comparison between the experiment and the model. 

5.2.1 Nominal atmospheric pressure 

The comparison between the predicted flow front assuming 101 kPa inlet pressure and the 
measured values (Panel 092904) is presented in Figures 58 and 59. For reference, two 
infusions of 1 stack of MAWK are also shown. The difference between the experiment and 
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the simulation is evident as 3DINFIL predicts an infiltration time of about 120 seconds 
compared to 240 seconds for the experiments. However, the flow front evolution in the 
medium predicted with 3DINFIL is in good agreement with a one stack MAWK preform, the 
same used for the skin of the hat stiffened panel (Figure 58). The agreement for the tool side 
between 3DINFIL and the one stack panel is not as good as the presence of the hat inserts 
affects the flow from the medium to the tool plate (Figure 59). It is reasonable to expect that 
the flow front at the tool interface of the hat stiffened panel be slower compared to a 1 stack 
preform. Therefore, the large difference between the 3DINFIL predictions and the actual 
experiment for the hat stiffened preform suggests that something was slowing down the 
infiltration. It could be possible that the resin line or the spiral tubing could be pinched 
creating a restriction to flow and consequently a pressure drop. 



Flow front position (m) 


Figure 58 Comparison between predicted and measured flow front position on the 

medium side. 
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Figure 59 Comparison between predicted and measured flow front position on the tool 

side. 

In order to investigate the effect of a pressure drop in the resin line on the preform infiltration 
behaviour, a series of simulation are conducted where the resin pressure is set to 80, 60 and 
50 kPa. The results are presented in Figure 60 and Figure 61. The effect of a decrease in the 
resin pressure also decreases the discrepancy between the simulations and the experiment. 
The 50 kPa resin pressure simulation agrees very well with the experiments. This tends to 
confirm that something in the resin line causes a significant pressure drop from atmospheric 
pressure. An indication of such a problem comes also from the very uneven flow front in the 
preform width observed in the early stages of the infiltration (Figure 54 frames 1 and 2). In 
typical 1 stack or even 4 stacks infiltration, the flow front is very unifonn across the preform 
width. The infusion time difference between the experiment and the model right under the 
hat stiffener inserts (Figure 57, Point 16) can be explained by the fact that the aluminium 
insert provides additional compaction on the preform as they act as pressure intensifies. 
This cannot be accounted for in the infusion simulation. 
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Distance (m) 


Figure 60 Predicted and measured flow front position at various pressures over media 

side. 



Distance (m) 


Figure 61 Predicted and measured flow front position at various pressures on tool side. 
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6 Conclusion 

Several conclusions can be drawn from this work and they are summarized below. 

• The VR56-19 resin have different cure kinetic and viscosity compared to the VR56-4 
system. 

• 3DINFIL can predict the infiltration of the uniaxial tackified preform, but great care 
should be taken when selecting the initial preform thickness as the infiltration time is 
very sensitive to this parameter particularly for this type of prefonn 

• A sensitivity analysis of the variation of material and processing conditions shows 
that the preform compaction behavior has a strong effect on the preform infiltration. 
The characterization of this material property is therefore very critical for accurate 
prediction of preform infiltration and compaction. 

• The infiltration of complex shapes is possible with 3DINFIL and the optimization of 
the infiltration point and medium length could control the part final fiber volume 
fraction. 

• The infiltration of a complex hat stiffened preform is simulated with 3DINFIL and 
results suggest that the flow during the actual experiment was affected by a restriction 
in the resin line leading to a significant pressure drop and un-uniform flow front 
across the preform width. 
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